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A spin-orbital superexchange Hamiltonian in a Mott insulator with t2g orbital degeneracy is inves- 
tigated. More specifically, we focus on a spin ferromagnetic state of the model and study a collective 
behavior of orbital angular momentum. Orbital order in the model occurs in a nontrivial way - 
it is stabilized exclusively by quantum effects through the order-from-disorder mechanism. Several 
energetically equivalent orbital orderings are identified. Some of them are specified by a quadrupole 
ordering and have no unquenched angular momentum at low energy. Other states correspond to 
a noncoUinear ordering of the orbital angular momentum and show the magnetic Bragg peaks at 
specific positions. Order parameters are unusually small because of strong quantum fiuctuations. 
Orbital contribution to the resonant x-ray scattering is discussed. The dynamical magnetic structure 
factor in difi'erent ordered states is calculated. Predictions made should help to observe elementary 
excitations of orbitals and also to identify the type of the orbital order in ferromagnetic titanates. 
Including further a relativistic spin-orbital coupling, we derive an effective low-energy spin Hamil- 
tonian and calculate a spin-wave spectrum, which is in good agreement with recent experimental 
observations in YTiOs. 



PACS numbers; 75.10.-b, 75.30.Ds, 75.30.Et, 78.70.Nx 



I. INTRODUCTION 

Many transition-metal oxides fall into the category of 
Mott insulators,^ in which the large degeneracy of atomic 
states remain unquenched down to low energies. Of par- 
ticular importance here is the role being played by orbital 
degeneracy inherent to perovskite lattices. An additional 
degeneracy of low-energy states and the extreme sensi- 
tivity of the chemical bonds to the spatial orientation of 
orbitals lead to frustrating interactions and a variety of 
competing phases that are tunable by moderate external 
fields.2 

As "orbital physics" has started to become an essential 
ingredient of the physics of transition-metal oxides, more 
efforts are necessary to develop quantum many-body the- 
ory of coupled spin-orbital systems in order to under- 
stand specific features of the orderings and fluctuations 
in these models. Earlier work has emphasized a "classical 
part" of the problem, focusing mainly on the strong in- 
terplay between classical spin and orbital configurations. 
It is implicitly assumed that at low temperature orbitals 
are frozen in a certain static pattern that optimizes both 
superexchange (SE) and orbital-lattice (JT) couphngs. 
Such a classical approach has been used with a great 
success as a theoretical guide in studies of magnetism of 
transition-metal oxides. 

Recent experimental developments indicate, however, 
the limitations of this standard picture. It has been ar- 
gued that quantum fluctuations of orbitals might some- 
times be of crucial importance, hence quantum version 
of the orbital physics is needed. New concepts, such as 
three-dimensional orbital liquid in LaTiOa (Refs. 3 and 
onlineciteKHAOO) and one-dimensional orbital chains 



showing Heisenberg-like orbital dynamics in cubic vana- 
dates (Ref. 5), have been proposed. It is not accidental 
that these ideas emerge from a study of titanates and 
vanadates having t\g and t\g electronic configurations, 
respectively. This is because of (i) large, threefold de- 
generacy, and (ii) a special rotational symmetry of t2g or- 
bitals. Another crucial point is that (iii) the JT coupling 
is relatively weak for t2g systems. Indeed, JT-like elon- 
gation of octahedra on titanates and vanadates is much 
smaller compared with typical JT distortions in mangan- 
ites with Cg orbitals. One may therefore think that t2g or- 
bital states are much less affected by electron-lattice cou- 
pling, and an intrinsic dynamics of coupled spin-orbital 
system governed by electronic superexchange interaction 
becomes the decisive factor in a first place. Effects of 
lattice distortions (which are always present) can then 
be accounted for in a next step. This point has actu- 
ally been emphasized long ago by Kugel and Khomskii^ 
indicating also very peculiar specific features of t2g spin- 
orbital models.^ 

The aim of this paper is twofold. First, we study the 
orbital state and orbital quantum dynamics in the ferro- 
magnetic state of superexchange model with t2g orbital 
degeneracy in a cubic lattice. Second, we discuss the re- 
sults in context of the magnetic properties of YTiOa, a 
rare example of a ferromagnetic Mott insulator.®'^ Recent 
spin-wave data shows that the ferromagnetic (F) state 
of this material is highly isotropic having the same ex- 
change couplings in all cubic directions.^" This is in sharp 
contrast with expectations from the conventional orbital 
ordering picture resulting commonly in a strong spatial 
anisotropy of the spin exchange bonds. ^^'^^ This obser- 
vation already indicates a rather unusual orbital state in 
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YTiOs- We would like also to understand a mechanism 
which stabilizes such a isotropic F-state in YTiOa, hav- 
ing in mind that its sister compound, LaTiOs shows a 
completely different, antiferromagnetic (AF) state. Cu- 
riously enough, spin-exchange couplings in LaTiOa are 
also of cubic symmetry and spin gap is also small, and 
these observations were understood in terms of fluctuat- 
ing orbitals.*'^^ 

Wc argue that AF and F states in t2g SE model arc 
actually very close in energy and strongly compete. This 
is because in both states there are large-scale orbital fluc- 
tuations gaining almost the same amount of the superex- 
change energy. Yet the AF state is slightly lower be- 
cause of an additional, composite spin-orbital fluctuation. 
However, an external parameter, namely, a larger distor- 
tion of Ti-O-Ti bonds due to a small size of Y-ion in case 
YTiOa induces an additional ferromagnetic coupling in 
all three directions, and stabilizes the spin F state. This 
distortion induces also a gap for orbital excitations. The 
orbital order pattern is very specific, and it supports ex- 
actly the same ferrocouplings in all three directions. We 
derive an effective spin Hamiltonian which includes ef- 
fects of the relativistic spin-orbital coupling as well, and 
show that this Hamiltonian leads to spin-wave dispersion 
and spin gap consistent with experimental observations. 
Some of these results were presented in Ref. 14. 

The following part of this paper is structured as fol- 
lows: Section H presents superexchange Hamiltonian in 
tig Mott insulator; Sec. HI presents orbital ordering and 
fluctuations in the F state; Sec. IV presents stabiliza- 
tion of the F state by Ti-O-Ti bond distortion; Sec. V 
presents orbital gap induced by Ti-O-Ti bond distor- 
tion; Effective spin Hamiltonian and magnons are given 
in Sec. (VI); Sec. VII lists predictions for resonant x- 
ray scattering; Sec. VIII gives orbital angular momen- 
tum contribution to the neutron scattering cross-section; 
(IX) Summary and discussion. The appendixes A-C con- 
tain some lengthy equations; Appendix D shows magnon 
softening by orbital fluctuations; Appendix E shows spin 
interactions in a previously reported orbital state for 
YTiOa. 



II. HAMILTONIAN 



tions leads to superexchange interactions, which in or- 
bitally degenerate systems strongly depends on the or- 
bital structure. In general, it can be written as 



H: 



SE 



Si-Si + -\J, 



?(7) 



2 *J 



(1) 



where the orbital operators J-^^ and K^^^ depend on 
bond directions 7(= a,b,c). In a t2g system like the 
titanates they are given by the following expressions:^'^ 



47^ = ^1 



i(n+r2)A(;)-i(r2-r3)B^ 



-jin - r2){ni+nj) 



(7) 



(7) 



i(ri-r2)A(7) + i(r2-r3)B, 



(7) 
ij 



-j{n +r2){ni + nj) 



(7) 



(2) 



(3) 



where Jse = ^t^/U. The coefficients ri = 1/(1-877), 
r2 = 1/(1 — 77), and r^, = 1/(1 + 2?;) originate from 
the Hund's splitting of the excited fl^g multiplet via 
ry = Jh/U. Reference 15 gives Jh ~ 0.64 eV and the 
multiplet averaged Coulomb interaction (= U — ^Jh) ~ 
4 eV, from which representative values U 5.4 eV and 
r] - 0.12 follow. 

The operators A^ij \ B[j \ and can conveniently be 
represented in terms of constrained particles (orbitons) 
Oi, hi, Ci with nia+nib+riic = 1 corresponding to t^g levels 
of yz, xz, xy symmetry, respectively. (This notation is 
motivated by the fact that each t2g orbital is orthogonal 
to one of the cubic axes a,b,c.) Namely, 

= niaUja + nibUjb + albib^jaj + bjajatftj, (4) 
B^f = niaUja + nibUjb + a\hia\hj + blaibjaj, (5) 
n^"^ = + riib (6) 



for the pair along the c axis. Similar expressions are ob- 
tained for the exchange bonds along the axes a and b, by 
replacing orbitons (a, 6) in Eqs. (4)- (6) by (6, c) and (c, a) 
pairs, respectively.^^ Another useful representation of or- 
bital exchange operators is via the angular momentum 
operators of t2g level,'' using the following relations: 



A. Superexchange interaction in orbital system 1^= ilc'b ~ b^ c) ,ly= i[a^ C — a) ,lz= i{b^ <l — a^b) . (7) 



Wc start with a discussion of the model Hamiltonian. 
In Mott insulators, the competition between kinetic and 
potential energies is resolved in favor of strong correla- 
tions that lead to a localized electron picture. Charge 
localization is however not perfect: electrons still make 
virtual excursions to neighboring sites in order to retain 
their kinetic energy at least partially. In terminology of 
Mott-Hubbard insulators,^ the zero-point charge motion 
is described as a high-energy virtual transition across the 
Mott gap. Kinetic energy associated with these transi- 



In terms of these angular momentum operators, a!.j\ 



b'i- \ and nP' are represented as 



.(7) 



A? = [(1 - ilW - + {iJyUiyQj] + [x-- y], (8) 

Blf - [(1 - llW - ll)j + {IJyUUyh] + [X^ yl (9) 



(10) 



Expressions of these operators for a and b bonds are given 
by replacing two component of the angular momentum 
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{lx,ly) in Eqs. (8)-(10) with (ly-lz) and (lz,lx), respec- 
tively. Angular and quadrupolc momentum representa- 
tion of the t2g superexchangc has recently been used also 
in Rcf. 17. In addition to Eqs. (4)-(5) and (8)-(9), it is 

also useful to represent A-J-* and b'^'J'' in terms of auxil- 
iary orbital pseudospins: 

4]^ = 2{n-^ + ^''\ (11) 

BS^=2{r,^r, + ^f. (12) 

Here t-''^ is a pseudospin one-half operating on the sub- 
space of orbital doublet {a, (3)'^'^^ active on a given 7- 
bond. Namely, pseudospin 7^'^^ operates on the subspace 

spanned by (a, b) pair of orbitons, while -7^^"'' and t^'''^ act 
on (6, c) and (c, a) doublets, respectively. A symbol (8) 
denotes a product Tj t^- = rf r? + {t^t^ + Ti^T~)/2. 

B. Ferromagnetic state 

The ferromagnetic state of a Mott insulator is usually 
thought due to a particular orbital ordering that opti- 
mizes the intraatomic Hund's exchange of electrons in 
doubly occupied virtual states. This is not the whole 
story, however. Neglect for a moment the Hund's cou- 
pling terms in Eqs. (2) and (3) (consider 77 ^ limit). 
The Hamiltonian obtains then the following structure: 

Ho = JseY.2 (Si ■ Sj + i) (ji ■ fj + ^mnj) . (13) 

(The unessential energy shift, —Jse, is not shown here). 
Regarding a single bond, one notices that spin coupling 
may be of cither sign, depending on the intersite orbital 
correlations. Singlet correlations of orbital pseudospins 
tend to align spins ferromagnetically, hence cooperating 
with Hund's rule effects. In systems with large, classical 
spins (e.g., vanadates), such a quantum orbital singlet 
controls the ground state. ^'^^ In quantum spin one-half 
case of titanates, however, (spin triplet) x (orbital singlet) 
and (spin singlet) x (orbital triplet) configurations are de- 
generate and compete. In a lattice, quantum resonances 
between these configurations are possible.^ In general, t2g 
superexchange Hamiltonian Eq. (1) represents a highly 
frustrated many-body problem. Wc will return to the 
interplay between antiferromagnetic and ferromagnetic 
states later on, while focusing now on the ferromagnetic 
state realized in YTiOa. 

In the spin saturated state, Eqs. (l)-(3) are simplified 
to: 

Horb = -riJsE + ^nJsE a'i]\ (14) 

{ij} 



where A\j is given by cither of Eqs. (4), (8), and (11). We 
consider the orbital order and dynamics in this Hamilto- 
nian. The effects of the dynamical coupling between spin 
excitations and orbitals that is present in Eq. (1) will also 
be discussed in the context of magnon spectra. 

III. ORBITAL ORDERING AND EXCITATIONS 
A. Discussion of possible orderings 

Even though spin as well as composite spin/orbital dy- 
namics is "switched off" in Hamiltonian Eq. (14), it still 
contains nontrivial physics. 

It is useful to look at the structure of Horb from dif- 
ferent points of view, (i) On a given bond, the operator 

A^J^ acts within a particular doublet of equivalent or- 
bitals. Spin-like physics, that is the formation of orbital 
singlets is therefore possible, (ii) On the other hand, 
interactions on different bonds are competing: they in- 
volve different doublets, thus frustrating each other. This 
brings about a Potts-model-like frustration, from which 
the high degeneracy of classical orbital configurations fol- 
lows, (iii) Finally, we observe in A^J"* a pseudospin I = 1 
interaction of pure biquadratic form [sec Eq. (8)]. Would 
I be a classical vector, it could change its sign at any 
site independently. Such a local (so-called Z2) symmcitry 
and the associated degeneracy of the classical states tell 
us that angular momentum ordering, if any, must be of 
pure quantum origin. 

The above points (i)-(iii) govern the underlying physics 
of the orbital Hamiltonian. We need to find such classi- 
cal states that provide best zero point energy when we 
switch on the quantum fluctuations. In other words, cer- 
tain classical orbital patterns will be selected and sta- 
bilized by quantum effects via the ordcr-from-disordcr 
mechanism. ^'^ Normally, these orderings are expected to 
be along symmetric orientations of the crystal depending 
on symmetry of the imder lying interactions. 

By inspection of the global structure of A^^^ [Eq. (8)], 
one observes that the non-cross terms, such as (1 — 
?^)i(l — lx)jy are definitely positive. However, the cross 
terms, {lxly)i{lylx)j and {lylx)i{lxly)j [which change the 
"color" of orbitals, see Eq. (4)], can be made negative on 
all the bonds simultaneously, if (i) on every bond, two par- 
ticular components of li and Ij are antiparallel, and (ii) 
remaining third components are parallel. For c bonds the 
rule reads as: kzljz and kxljx are both negative, while ly 
components are parallel. (In terms of orbitons: Ci and cj 
are in antiphase, and well; but 6j and bj have the 
same phase.) We find only two topologically different ar- 
rangements [called (a) and (b)], which can accommodate 
this curious mixture of "2/3 antifcrro" plus "1/3 ferro" 
correlations (see Figs. 1 and 2). In the state (a), sub- 
lattice unit vectors are along the cubic diagonals [111], 
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(a) (b) 

FIG. 1. Two kinds of four-sublatticc structure for orbital 
orderines in SDin ferrornaEuotic tin suDerexchange model. 
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(a) (b) 

FIG. 2. Arrangement of the local quantization axes in 
states (a) and (b). Arrows indicate the quantization axes 
at each site, and represent also a snapshot of local correlar- 
tions of angular momentum: on every bond, two out of three 
components of I are correlated antiparallel. 

while in state (b), sublattice unit vectors are [110] and 
[002]. 

For technical reasons, it is useful to introduce new 

quantization axes. This is done in two steps. First, we in- 
troduce local, sublattice specified quantization axes (see 
Fig. 2): 



ix,y,z) 

{x,y,z) 
(x,y,z) 
(x, y, z) 



{x,y,z), 
{-x,-y,z), 

{-x,y,-z), 



(15) 
lia and or- 



After corresponding sign transformations of 
bitons, one obtains 

A\f = niaUja + rhibTijb - albib^jaj - blaiO^jbj (16) 

= [(1 - llW - iDo - iUyUlyQj] + y]. (17) 

^From now on, a sublattice structure will not enter in 
the excitation spectrum. From the above observations 
it is also clear that all the components of I are equally 
needed to optimize all the three directions. We antic- 
ipate therefore that the cubic diagonals are "easy" (or 
"hard") axes for I fluctuations/orderings (recall that the 
Hamiltonian has no rotational symmetry for I vector). 
Therefore, it is convenient to further rotate the quanti- 
zation axis so that new z-axis (denoted as z) corresponds 
to [111] direction. This is done as follows: 



ll — R 



where U = {kx^ky, kz) and R is given by 



1 



s c — s 1 



R = —= c — s c + s 1 
^\ -1 -11 



(18) 



(19) 



with c = 1/2 and s = Here, new x- and y-axes are 

taken to be symmetric with respect to the [110] direction. 
Annihilation operators for constrained particles obey the 
same transformation. 



(20) 



Explicit expressions for the wave functions Va are ob- 
tained by reversing Eq. (20) as follows: 




-^[{C + S)dyz + (C - S)rfa: 

-^[(c - s)dy^ + (c + s)dx 



d 



xv\i 



(21) 



In Fig. 3, we show schematic pictures of these orbitals. 

By construction, tfj^ is symmetric with respect to the ro- 
tation around [111] direction, having simply iz? — r^ sym- 
metry, while ij^a and -i/ij^ are symmetric with respect to 
[110] axis. At the end, the orbital Hamiltonian in a ro- 
tated basis obtains the following form (symbol "tilde" de- 
noting rotated axes is implied for angular and quadrupole 
operators below, and constant energy shift is dropped 
out): 



{ij) 



(22) 



with 



(7) 



2 1 

— (1 QizQjz) 2^iz^jz 

+ 2 + T-iQx — TqTi — TiTo)ij 



-{QxQx-ToTo + T_iT_,)\]^ 
2^ ,„ „ (^) 2, 



3 

'^7;hzilx + ly)j 



2 

'niJ'x iy)i \^x~^iy)j 



2 i^x + iy)i^^h^ 



(23) 



where Qz and represent the quadrupole moment op- 



erators with eg symmetry, 3z 



and 



y , respec- 



tively. To = and r±i = Ty ± T^, where T^, Ty, and 



4 



FIG. 3. (Color online). Orbitals in a new basis xyz speci- 
fied by transformation Eq. (19). 



Tx represent the quadrupole moment operators with t2g 
symmetry of xy, xz^ and yz, respectively. These oper- 
ators are expressed in terms of angular momentum and 
orbiton operators as follows: 







Uc - 


1. 




-l^ -l^ ~ 

~ '■x '•y ~ "b 


■Ha, 






— ^x^y ^y^x — 




+ 6^ a), 




^y^z ^ ^z^y ^ 






Ty 


— ^x^z ^z^x — 


-(2ta 





(24) 

7-dependence of quadrupole moment operators is ob- 
tained by changing l^.y in Eq. (24) to 

^(7) _ f -cl^ ± sly, for 7 = a(6) 
" ^ Ix, for 7 = c 



7(7) _ 



-clyTslx, for 7 = a(6) 



for 7 



(25) 



Explicit expressions for Q^a'^ and T^'''-* are given in Ap- 
pendix A. It should be noted that, among eight opera- 
tors, namely, five quadrupole moment and three angular 
momentum operators, only four operators are indepen- 
dent of each other because of the local constraint among 
orbiton operators. 

Although it looks a bit complicated, the rotated Hamil- 
tonian obtains a well-structured form. The first and sec- 
ond terms of Eq. (23) represent "Ising"-like interaction 
for quadrupole moments and angular momenta. This 
part of the Hamiltonian stabilizes the ordering (conden- 
sation) of an appropriate orbiton. On the other hand, 
the other terms represent fluctuations of Qx, T^, Ty, 
and transverse components of angular momenta and 
ly. These terms generate dispersion of the orbital exci- 
tations. 

Ordered states, promoted by the "Ising" part of inter- 
actions, can be characterized by the quadrupole moment 



I n 

FIG. 4. (Color online). Schematic pictures of the orbital 
orderings. Left: Real orbital ordering 1(a). Right: Complex 
orbital ordering 11(a). Here, the absolute values of the wave 
functions are presented. Arrows represent the directions of 
angular magnetic momenta in the orbital magnetic state 11(a). 



Q ~ {Qz) (Q-order may couple to a lattice distortion 
of _D3(; symmetry), and the angular magnetic moment 
™; = {^z}- We notice that the magnetic, hzljz term in 
the first line of Eq. (23) is generated by quantum commu- 
tation rules when we rotate Horb\ this makes explicit that 
the Z2 symmetry is only a classical one and emphasizes 
the quantum origin of orbital magnetism. 

As it follows from the definition of Qz, quadrupole or- 
dering with finite Q but zero 1^ corresponds to a con- 
densation of the c orbiton. We call this solution state I. 
Classically, Q = 1 in this state. On the other hand, con- 
densation of the complex orbital (d — ib) / \/2 generates 
a finite magnetic moment (to/ = 1 classically), the state 
called II. The orbital patterns in a classical states I and 
II are shown in Fig. 4. We now focus on fluctuations of 
orbitals, and show that the excitation spectra are in fact 
identical in these states. Moreover, we will obtain that 
the states I and II can smoothly be connected by a con- 
tinuous phase rotation of the condensate wave function. 
Noticing that an arbitrary cubic diagonal could be taken 
as z and having in mind also two structures in Fig. 1, one 
obtains a multitude of degenerate states. This makes, in 
fact, all the orderings very fragile. 



B. Orbital quadrupole order 

This state is driven by a condensation of c orbital, that 
is, ordering of the orbital t/)^ in Eq.(21). To obtain a 
linear orbital wave Hamiltonian, we resolve a constraint 
^1 - na 



as c — c) — 



I — ni, and expand Eq. (23) up to 
second order in d and b. The result is (in units oiriJss)'- 

How = ^{nm + 

i 

where z = 6, and 
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(000) 
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FIG. 5. Upper panel: Orbiton dispersions (in units of 
riJss), obtained in a linear spin- wave approximation. Lower 
panel: Positions of soft modes are shown by the thick lines. 



T±i = -(at+a)T(bt + 6), 
lx = i{b — b^), ly = i{a^ — a). 



(27) 



In a momentum space, this linearized Hamiltonian reads 
as 

How = J2 Kfc + + ^(71 + 72)(aia^_g + a^a_^) 

k 

+^(7i-72)(6p_^^ + 6,^6_g) 

-73(46Lj + %4_g)], (28) 

where 71, 72, and 73 are defined as 71 = {cx + Cy + Cz)/3, 
12 = \/3(cy - Cx)/6, and 73 = (20^ - - Cy)/6, respec- 
tively, with c„ = cos ka ■ The Hamiltonian is diagonalized 
by using Bogoliubov transformation (see for details Ap- 
pendix B). One obtains: 

How = Yl ('^Ife + '^2fc"2fc°=2fc) + ^0 , (29) 



where w^g = ^1 - (71 -|- k)^ and lo^j^ = y'l - (71 - 

with K = 1/7! + 73. The dispersion relations of orbitons 

arc presented in Fig. 5. Orbital excitations arc char- 
acterized by the flat dispersion with zero energy along 
{0,0, kz), {TT,Tr,kz) and their equivalent directions. As 
discussed later, interaction effects open the gap along 
(tt, TT, ^2) and equivalent ones, while zero modes along 
{kx,0,0), {0,ky,0), and {0,0, kz) are protected by the 
underlying symmetry of the model. 

A constant Eq in Eq. (29) represents the energy gain 
due to the quantum fluctuations. It is given by 



^0 = 1 + <^2k) - 1 = -0-214 {nJsE). 



(30) 



We may compare this result with ground-state en- 
ergy of the orbital disordered AF state: Eq = 
— 0.33^^^^^^ J SE,"^ where the result of Ref. 4 is corrected 
for the finite values of rj. For realistic values of the Hund's 
couphng, say rj = Jh/U = 0.12, this gives Eq = —0.285 
(in units of ViJse) in AF state. It is noticed that fer- 
romagnetic and AF states are almost degenerate. Still, 
the ferromagnetic state is higher than the AF state, so 
its stabilization in YTi03 requires an additional effects 
as discussed in Sec. IV. 

Due to the flat mode, one may expect strong orbital 
fluctuations in the ground state. Indeed, number of the 
excited bosons d and b is large even at T = 0: 

(n,a+n,^)=-l + jEf^ + ;Z:) =°-^^- ^^^^ 



ife 



UJ. 



2k 



This reduces the condensate density to {riic) = 0.46. 
Consequently, the quadrupole order parameter is ob- 
tained to be rather small: Q = 0.19. Reduction of 
quadrupole order Q implies that electron density is much 
less anisotropic than that shown for the classical state in 
Fig. 4(a). Including fluctuation effects, that is, finite pop- 
ulation of d and b orbitals, electron density at site 1 is 
given by pi{r) = ngV'i + natpl + niip?. Using Eq. (21), 
one then finds 



di 



^xy) 



+ -^Q{dyzdx 



dyzdxy "4~ dxzdxy)' 



(32) 



dyzdxy 



'dyzdxz ~^ dyzdxy dxzdxy^ for site 3, 



dxzdxy) for site 4. In Fig. 6(b), 



Electron density at other sites is given by a similar equa- 
tion, where the second, of t2g symmetry term, is different 
for different sublattices. Namely, it is {dyzdxz — dyzdxy — 
dxzdxy) for site 2, (- 
and {-dyzdxz 

we present the electron distribution given by Eq. (32). 
For comparison, we show in Fig. 6(a) the electron distri- 
bution where a, b, and c are equally occupied. At finite 
Q, the electron density pi{f) is slightly elongated along 
[111] direction. Thus, we expect the quadrupole ordered 
state to be further stabilized by the electron lattice cou- 
pling, although this coupling is expected to be weak for 
t2g orbitals. 

The anomalous reduction of the order parameter is 
due to the highly frustrated nature of the interactions 
in Eq. (14). A special, non- spin- like feature of all orbital 
models is that orbitals are bond selective, resulting in a 
pathological degeneracy of classical states. This leads to 
soft modes [observe that is just flat along (0, 0, tt) 

and equivalent directions]. These soft modes have their 
origin in special symmetry properties of the t2g orbital 
model Eq. (14), which result in conservation laws with 
important consequences. Namely, the total number of 
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(a) (b) (c) 

FIG. 6. (Color online). Spatial electron distribut: 
different states, (a) Disordered state, = «(, = 

(b) Quadrupole ordered state, ng = 0.46 and = Wj; = 
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orbitals of each "color" (a, b, c) are conserved during su- 
perexchange process, as can easily be seen from Eq.(4). 
Moreover, as i2g-orbitals can hop only along two direc- 
tions [say, xy-orbital motion is restricted to (ab) planes] , 
the orbital number is conserved on each plane separately. 
Formally, these conservation rules are reflected by a pos- 
sibility of uniform phase transformation of orbiton op- 
erators, e.g., a — > a exp{i<f>a), etc., which leaves the or- 
bital Hamiltonian invariant. These continuous symme- 
tries are spontaneously broken in above ordered states. 
The breaking of continuous symmetry is usually followed 
by the generation of gapless Goldstone modes. This is 
precisely what happens in the t2g orbitally degenerate 
model. In fact, soft modes obtained above have two- 
dimensional (2D) feature (stemming from 2D geometry 
of the t2g orbital hoppings). As a result, long range or- 
bital order is possible only at zero temperature, just like 
in 2D Heisenberg models. Formally, this is manifested 
as a divergence (as In y ) of the number of thermally ex- 
cited orbitons, {ria + ni), if one attempts to calculate this 
quantity at finite temperature, including the Bosc popu- 
lation factor in Eq.(31). Soft modes will be discussed in 
more detail later on. 



C. Orbital magnetic order 

In order to describe the magnetic ordering (denoted 
above by state II), let us introduce orbital states 



■^{d-ib), b^^{h + ib). 



(33) 



In these states, angular momentum has eigenvalues Iz = 
±1 (on local axes), respectively. A condensation of a 
leads therefore to the magnetic ordering (with Iz = I clas- 
sically), and it is associated with ordering of the complex 
orbital: 



e*3 +da:ze '3-d 



(34) 



On global axes, this order is noncoUincar, as shown in 
Fig. 4(b). We consider fluctuations of this state. Using 



the relation a = = y/1 — ni ~ ric and expanding in- 
teractions in Eq. (23) up to second order in b and c, we 
obtain the following linearized Hamiltonian: 

How = ^{na + riit) 

i 

+ 2^ I^Qix?} _i + - iQjx — TioTii — TiiTjo 



~^Qix Qjx 



rp(l) rp(l) 



(35) 



Here, the operators Q, and T arc linear functions of b 
and c: 



+ C) , ly 



T. 



±1 



(5t + 6), To = -i(6t-6), 



(36) 



We introduce now new operators a, (3: 



a = ^{b + e^n), /3 = -L(6-e''^c), (37) 

with (y5 = 7r/4. Remarkably, after this transformation 
the linearized Hamiltonian obtains in a momentum space 
the same form as Eq. (28), where a, b are just replaced by 
a, (3, and 72 and 73 are interchanged. Exchange of 72 and 
73 does not affect the excitation spectrum, however, as 
they enter in oj^ ^j: via the parameter k = vtI+Ts only. 
Thus, the ground-state energy in the orbital magnetic 
state Eq is given again by Eq. (30), so the states I and 
II are degenerate even on quantum level. Similarly, the 
number of out of condensate bosons is also obtained from 
Eq. (31). This gives {ua + np) = 0.54, and values for 
the angular and quadrupole momentum order parameters 
follow: 



mi = l- -{rio 



np) =0.19, 

+ nf,) = -0.095. 



(38) 



Electron-density distribution in the state II, given by 
Eq. (32) with above value of Q, is shown in Fig. 6 (c). 
Electron cloud in the magnetic state is almost of cubic 
symmetry, being just slightly contracted along [111] di- 
rection (opposite to that in the state I). Thus the energy 
gain from the orbital lattice coupling is smaller in this 
state. In principle, orbital magnetic order could be sup- 
ported by a relativistic spin-orbital coupling; however, in 
the noncoUinear state driven by the superexchange inter- 
action, uniform component of the orbital moment is zero 
[see Fig. 4(b)], hence the coupling to the spin ferromag- 
netism vanishes in linear order. The spin-orbital energy 
gain in a second order is possible though, via the canting 
of spins towards orbital magnetic pattern. 



7 



D. Soft modes: density-phase formulation 

Having obtained an identical excitation spectrum for 
states with apparently different ground-state conden- 
sates, we would like to unify these states. It is conve- 
nient to use a different approach, that is, the density- 
phase formulation by Popov (Ref. 20), nowadays called 
"radial gauge" . This formalism is particularly useful also 
to clarify the physical origin of the soft modes obtained 
above. In the radial gauge, density and phase degrees of 
freedom of the constrained particles are emphasized. We 
work in a basis obtained by first transformation, Eq.(15), 
and represent the orbiton operators entering in Eq.(16) 
as follows: 



/piaC*^'" (a = a, b, c). 



(39) 



Further, the density and phase degrees of freedom are 
parametrized as 



and 



Pia{b) = Po + -^{cri±sXi), 



Po 



ia{b) =^i+ Cipi ± sQi 
die ~ — fi; 



(40) 



(41) 



respectively. Here, po{= 1/3) is an average electron den- 
sity on each orbital. The phase Oj, common to all the 
three orbitons, can as usually be absorbed by the con- 
straint field, while the local constraint itself is explicitly 
resolved by parametrization (40). The physically active 
degrees of freedom are therefore r and A fields for the 
amplitude fluctuations, and 9 and ip for the phase fluctu- 
ations. We recall that the coefficients c = 1/2, s = V3/2. 

To start with, let us neglect for a moment the am- 
plitude fluctuations, and focus on the phase dependence 
of the classical condcinsate wave function. In terms of 
phases in Eq.(41), it is written as follows (up to unessen- 
tial overall phase factor): 

V(e,<p) = V^{(i,.e^(3^^+^') + d,,e^(3«,-.e) ^ ^^^Y (42) 

Here, we suppressed site dependence of the phases ip, 0, 
discarding for a while slow space variations of the conden- 
sate. Now, it is noticed that the quadrupole and mag- 
netic ordcrings [see Eqs.(21) and (34)] do follow from 
Eq.(42) when ip = = 0, and <^ = ^, = gj, respec- 
tively. Next observation is the orbital "color" conserva- 
tion rule in the ferromagnetic state. In the radial gauge, 
it is evident from Eq.(16), that the interactions do de- 
pend on the difference of the orbiton phases only, that is 
on Oia — 0ja, etc., so we can uniformly rotate the conden- 
sate function (42) by arbitrary phases ip, 6 with no energy 
cost. By such rotations, we can in fact mix quadrupole 



(state I) and magnetic (state II) orderings. Slow phase 
rotation of the condensate is precisely the origin of the 
soft modes obtained above. Because of the two dimen- 
sionality of t2g orbitals, the phases can spontaneously be 
fixed at zero temperature only. Of course, orbital-lattice 
and/or spin-orbital couplings may fix the phases, thus 
selecting a particular state even at finite T. 

We now turn to the excitations of the model in density- 
phase formulation. A^J^ in Eq. (16) is expressed as 



{y^PiaPja - ^PibP]b) 



+ '^y/PiaPjay/PibPjbi^ - C0S( 



<^f)}, (43) 



where 



— dia — 



^Ij'' for 7 = a and b bonds are 



given by replacing (a, b) in Eq. (43) with (6, c) and (c, a), 
respectively. In terms of the relevant phase degrees of 
freedom ip and 9, we obtain 



^ f \/3(-c^j ± s<^i), for 7 = a{b) 
\ ^/^ei. for 7 = c 



(44) 



The density operators pia are functionals of the and Ai 
fields. We may expand now the operator A^^^ in terms 
of the amplitude and phase variables r, A, 9. Keeping 
quadratic only terms in the expansion, one arrives at the 
following linearized Lagrangian for the phase-amplitude 
fluctuations 



Lorb — 



)+He^ + Hxr, (45) 



where uj is the Matsubara frequency. The first term in 
this equation originat(^s from the time derivative, kine- 
matic part of the Lagrangian —Ylia'-''l'3r'^i' P'^'^" 
duces dynamical coupling between the density and phase 
variables. Hg^ represents the phase fiuctuations, and it 
is obtained from the expansion of the second term in 
Eq. (43): 



L 2 



(ij) 



a(7) 



(46) 



The fluctuations of the condensate density about po are 
controlled by the Hamiltonian H\r, which in a linear ap- 
proximation reads as 
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,2 'ST 



H\r = -Po ^ 

{ij) 



(47) 



Here, 6^^^ is the differ ence between the electron den- 
sities on orbitals that are active on a given direction: 
(5 = Pj-^^ — p^^K Using parametrization (40), they are 
expressed via r and A fields as follows: 



■^{-cXi± sn), for7 = a(6) 



-^Xi. for 7 = c 



(48) 
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In a momentum space, Hq^ and H\r are represented 
as follows: 

He^ = + \<^M'' - 72¥'« (49) 

= II{^aA|A,H' + ia^lrg-l^ +72A,-r_,-}, (50) 

where = 1 - (71 + 73), a<p = 1 - (71 - 73), a\ = 
1 + (71 + 73), and = 1 + (71 - 73). 

Equations (45) and (49), (50) determine the orbital dy- 
namics in a harmonic approximation, which is equivalent 
to the previous linear orbital wave approach. Indeed, the 
obtained quadratic form in Eq.(45) can easily be diago- 
nalizcd giving exactly the same excitation spectrum, that 
is — (71 ± k)^ found in the preceding section. In ad- 
dition, the origin of zero energy excitations can clearly 
be identified now. In a classical limit [neglect dynami- 
cal term in Eq.(45)], the quadratic forms Hg^ (49) and 
H\r (50) can be diagonalized separately resulting in nor- 
mal modes with energies w^{k) ~ (1 — 71 ± k) in phase 
sector, and (fc) = (1 + 71 ± k) for the amplitude vari- 
ables, ivf vanishes on lines (kx, 0, 0), (0, ky, 0), (0, 0, kz) 
(sec Fig. 5). Therefore, zero-energy excitations on these 
lines do correspond to uniform phase rotations of or- 
bitons on different planes as discussed before. On the 
other hand, the normal modes Lo^{k) that are associated 
with the density of orbital occupancies possess zero lines 
at {kx, TT, tt), (tt, fcy , tt), (tt, tt, kz)- This reflects softness of 
the staggered fluctuations of orbitals [notice also that a 

uniform component, that is, S'l^'^ = S^'^^ +Sj'^\ only enters 
in Eq.(47)]. However, such soft modes in the amplitude 
sector arc not protected by physical conservation rules; 
therefore, they are expected to acquire a finite mass due 
to interaction effects that go beyond linear orbital wave 
approximation. To see this, one should consider unhar- 
monic terms in the expansion of Eq.(43). Most relevant 
term in that expansion is the interaction between the 
staggered fluctuations, that is, si''^ = sl'*^ — 6j^\ with 

uniform components 5u^. When such a term is kept, 
Eq.(47) is replaced by 

^lplT.{(^ + '^)^u+euS',f\ (51) 

(ij) 

The Hartree decoupling is applied here to the interac- 
tion term, with Ss,u = (l/4po)^('^s u)- These expectation 
values are finite due to the presence of quantum fluctua- 
tions in the ground state that are particularly enhanced 
in a staggered channel. In a momentum space, the above 
equation reads as 

Hxr = Z, ^{laA|Ag-|2 + ia^|r,-|2 + 7IA,- r_,-}, (52) 



where ax = 1 + (7i +7f ) and dr = l + {jf- 7f ). 7^ (n = 
1,2,3) are defined as jf^ = (1 — 2e)7„ with e — Su/Z^, 
and Zg; = (I + Eu + Ss) is an overall rescaling factor. 

Using now Hxr given by Eq.(52) in Lgrb [Eq. (45)], one 
obtains finally the following two eigenfrequencies: 

a;±(fc) = {1 - (1 - 2e)(7i ± - 2s{^, ± k)}'^' (53) 

in units of ^/Z^ r\JsE- It should be noted that a;± re- 
cover the orbiton energy 101^2 when e = 0. Using bare 
orbiton dispersions, the Hartree decoupling parameters 
es,u are calculated as follows: 

X i , (54) 

Numerical calculation gives Eg = 1-72 and Su = 0.59, 
reflecting that staggered fluctuations of densities are 
stronger. Thus we obtain ^fZl = 1.82 and e ~ 0.18. 
Dispersion relations of the orbital excitations, Eq.(53) 
calculated using these parameters are presented in Fig. 7 
(in units of ^/Z^ riJsE)- Staggered density fluctuations 
are gapped, and we are left now with true Goldstone 
phase modes protected by the symmetry of interactions. 



IV. WHY YTiOa HAS A FERROMAGNETIC 
GROUND STATE 

So far, we discussed t2g orbital physics on an ideal cu- 
bic lattice assuming a spin saturated ferromagnetic state. 
In the remainder of the paper, we apply the theory to 
the ferromagnetic state of Mott insulator YTi03. This 
requires some modiflcations of the theory implementing 
a specific feature of this material. On empirical grounds, 
it is well documented that Ti-O-Ti bond angle is an im- 
portant parameter controlling magnetic properties of ti- 
tanates i?Ti03.* The bond angle 9 gradually decreases 
from ~ 157 deg in LaTiOs to ~ 142 deg in YTiOa, due to 
lanthanum contraction effect that results in deviations of 
the lattice from an ideal perovskite structure. It is quite 
remarkable that such a small variation of the bond angle, 
driven by _R- ionic size effect, affects the magnetic state 
dramatically: It changes from isotropic AF as observed 
in LaTiOa to the isotropic ferromagnetic state in a Y- 
based compound, indicating strong competition between 
AF and F interactions in titanates. 

Ti-O-Ti bond distortion is important because it in- 
duces an unfrustrated ferromagnetic interaction, chang- 
ing thereby a delicate balance between AF and Ferro cou- 
plings that dynamically coexist and compete in ideal t2g 
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supcrcxchangc models like in Eq.(13). It was found in 
Sec. IIIB that the ferromagnetic state is slightly higher 
in energy that G-type AF one; the situation is however 
reversed when the bond angle is reduced below some crit- 
ical value, as we argue below. 

The bond distortion brings about the following two ef- 
fects, (i) Reduction of transfer intensity between nearst- 
neighbor (NN) t2g orbitals as t = A^^t^^^cos^ = 
toCOsO.^ Here, tdp-n{tdpa) is the transfer between Ti id 
and O 2p orbitals on the 7r(cr) bond, and A^p is the level 
difference between Ti 3d and O 2p states. Superexchange 



energy scale is then reduced as Jse = Jg^cos^^ with 
Jg!^ = At^/U . (ii) Generation of transfer intensity be- 
tween NN t2g and Cg orbitals, t' = A'^^tdpatdp-K sin 0. 
This transfer induces an additional SE interaction. For 
c-bond, we consider the following transfer term: 



H 



/(c) 



t' 



+ h.c.) + {i ^ j) 



(56) 



denotes the creation operator of electron in the eg or- 
bital with — symmetry. H' t for o- and 6-bonds are 
given by replacing {a, a) in Eq. (56) by [— (a— \/3/3)/2, 6] 
and [—(a + •\/3/3)/2, c], respectively, where (3 denotes 
the electron annihilation operator in the Cg orbital with 
— symmetry. By the second-order perturbation with 
respect to ff^, one obtains the SE interaction between NN 
t2g electrons. Here, energy of the intermediate cP excited 
states with spin triplet and singlet states between Cg and 
t2g electrons is given by U — 'iJH + ^cr and U— Jn + ^cr, 
respectively, where Acr is a cubic crystal-field splitting 



between Cg and t2g levels (so-called lODq). Explicit ex- 
pression for the new SE interaction is 



H 



/(c) 
SE 



'U 



1 



C/(l-377)(l-7?) 

2-3fi + AfiSi ■ Sj) {ilia + rlja), 



(57) 



where U = U + Ac,- and fj = Jh/U. H'se for a{b) 
bonds are given by replacing Ua in Eq. (57) by Uh (ric). 
It is stressed that this SE interaction is of the ferromag- 
netic sign, because t2g and eg orbitals are of the different 
symmetry, and the Hund coupling between them favors 
spin triplet state. As Ti-O-Ti bond angles in a, 6, and c- 
directions are almost the same,^ this interaction supports 
ferromagnetism equally in all three directions. Here we 
differ from Ref. 21, which considers t2g — eg hopping chan- 
nel along the c-axis only. 

Either in the orbital-ordered states [I (a) and I (b)] 
with f/ig = -^{dyz+dxz+dxy), or in the orbital-liquid one, 
average occupation number of each t2g orbital is given by 
Ua = rih = ric = 1/3. Thus, the spin interaction (57) in 



these orbital states becomes H'^pin 
with 



(58) 



One should notice that J' is proportional to sin^ 9 (via 
t'), and contains also the small number i] = Jh/U. This 
is because J' is caused by the Hund coupling between t2g 
and eg electrons in the virtually excited d'^{e^gt\g) state, 
which is evoked in SE process only in the presence of 
Ti-O-Ti bond angle distortion. 

The energy difference between AF and ferromagnetic 
phases stemming from is given by AE'g^ = ij', 

while that from Hse is given by AEse = (-0.33^^i|^ + 
0.2Uri)JsE (Sec. IIIB). The total SE energy difference 
between AF and ferromagnetic phases is then estimated 
as 



AE = 



-0.33 



n + r2 



1 ft'V/U^ ^ 



Jse. 



(59) 



Thus, AE obtains the following 9 dependence (at repre- 
sentative value T] = 0.12 for Hund's coupling parameter): 



AE 



-0.111 cos^ 9+U^^y( Z)\sm^ 9 



2 \tdp7r 



u 



jiO) 
'SE- 



(60) 



With the realistic parameters (typically, Acr is about 2 
eV) U/Acr = 2.5 and tdpa/tdp-K = 2, the transition from 
AF orbital-liquid phase to ferromagnetic orbital-ordered 
one occurs at the critical angle 6^'' = 136 deg in the 
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FIG. 8. (a) t2g-eg transfer originating from Ti-O-Ti 
bond distortion. Ti-O-Ti bond angle is denoted by 
d. (b) Energy difference AE (in units of J^^) between 
spin-AF(G)/orbital-liquid and spin-F/orbital-ordered states 
as function of 9. The solid line is a result of purely electronic, 
SE interactions. The broken line includes small JT energy 
gain SEjT in the orbitally ordered ferromagnetic state. Pa- 
rameters are ri = 0.12, U/Acr = 2.5, tdpcr/tdp-n = 2, and 
SEjT / ~ —0.04. (c) Schematic energy diagrams sd, 9 > 8c 
(left) and at 9 < 9c (right). 



SE model. This angle is slightly smaller than that ob- 
served in YTiOa {9 ^ 142 deg). Further, orbital-ordered 
state should be favored over the AF orbital-liquid state 
by orbital-lattice coupling. We simulate this by adding 
JT energy gain SEjt{< 0) to the energy of the ferro- 
magnetic orbital-ordered state. As shown in Fig. 8 (b), 
this increases 9c- A value of 6Ejt, which is required to 
obtain a realistic value 6c — 146 deg for titanates is small 
(—0.047^*^), so it might be hard to observe the associated 
^3(i-type distortion. 



V. ORBITAL GAP IN YTIO3 



Effect of TiOe tilting: In addition to the finite transfer 
between NN t2g-eg orbitals, octahedron tilting changes 
also the symmetry of NN t2g-t2g hopping matrix, mak- 
ing possible finite electron transfer between the NN or- 
bitals with different symmetry.^^ We show now that such 
a hopping leads to an important modification of the or- 



bital excitation spectrum, removing gapless Goldstone 
modes. Taking into account nondiagonal hoppings be- 
tween orbitals active on a given direction 



t" 



{alpj + f3la,) + {i ^ j) 



(7) 



(61) 



up to second order in small ratio t" /t, we obtain the fol- 
lowing correction to the interaction between NN orbitals 
in a spin-ferromagnetic state: 



H 



11(1) 

SE 



ri J; 



SE 



.h) 



(1 



(1) (7)1 

n- 'ny'\ 



(62) 



Here is given by either of Eqs. (5), (9), and (12). 

A crucial point is that this operator violates the orbital 
"color" conservation rule even in the fully spin polarized 
state. Therefore, uniform phase rotations (separately on 
each orbital fiavor) are not longer possible, hence the 
relative phases will be fixed and orbital gap will be gen- 
erated. 

Let us apply a radial gauge description, and focus on 
phase fluctuations, as the amplitude fluctuations have 
a large gap anyhow. In a local coordinates [defined by 
Eq.(15)], we find 



H 



nil) 
dip 



1 



riJsEpl 



t" 



E 



a(7) 



(63) 



which 



(7) 



Here, we expanded the function cos( 

enters in -B-^\ about 9,if — (this corresponds to the 

quadrupole ordered state, which is in fact favored by B 
term) . 

Effect of trigonal TiOe distortion: Our orbital state 
would be supported also by a trigonal (Dsd) distortion of 
TiOg. When this distortion is treated as a static one, an 
orbital feels the following potential: Hjt = —"^lEjrlQiz, 
where Ejt < represents the total JT energy gain. In a 
radial gauge, we obtain the following phase Hamiltonian 



(64) 



Summing up all contributions, Hg^ [Eq.(46)], H"g^, 
and Hg^ , one obtains the following phase Hamiltonian: 



Tjtot TT I TtH 



JT 



6ip 



(65) 



where Z 



y = (1 + (*^)2 + 'i\EjT\) is an overall factor. 

Here ag — \ — {j( + 73) and = 1 — (7/ — 73 ), with 
modified form factors 7^ = (1 — 2/)7„ (n = 1, 2, 3) where 

' t" 



f = {C-y 



\EjT\}/Zf. 
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Using now Eqs. (65) and (52) for Hg^p and H\r, re- 
spectively, we obtciin from Lorb 

[Eq. (45)] the following 

excitation spectrum: 

uj±{k) = {l - (1 - 2£)(1 - 2/)(7i ± «)2 

1 1/2 

-2(e-/)(7i±Ac)| . (66) 
This is given in units of Wort, which is defined as follows: 

Worb - V^IZj nJsE- (67) 

Worb represents the overall energy scale for orbital 
fluctuations in the problem. E xcitation gaps at 
(0,0,0) and (7r,7r, tt) are given by 2-\//(l — e) Worb and 
2y^£(l - f)Worb, respectively. Taking t"/t = 0.2 and 
\Ejt\ = 0.04(ri Js_e) as a representative values, we ob- 
tain / = 0.086. With e = 0.18 and y/Z^ = 1.82 estimated 
above (Sec. HID), one obtains Worb — 1-96 tiJse, and 
the lowest gap about 0.53Worb is then expected at (0, 0, 0) 
point. Thus, we expect that the orbital excitations in the 
modified model for YTiOs cover the energy window from 
~ nJsE to ~ 2 riJsE- 

VI. EFFECTIVE SPIN HAMILTONIAN 

The spin wave spectrum in YTiOs shows the "cubic 
symmetry" of the Heisenbcrg spin couplings: Ja — Jb — 
Jc.^° The magnon gap was found to be very small, al- 
most two orders of magnitude smaller that the magnon 
bandwidth (^ 20 mcV). It has been noticed that such 
an apparent simplicity of spin excitations, showing high 
isotropy in both real and spin spaces, is remarkable 
and puts strong constraints on possible orbital orderings. 
Spin wave excitations are examined in this section. Be- 
ing a test case for the above theory for orbitals in YTiOs, 
a comparison with experiment gives also an opportunity 
to estimate SE energy scale Jse in the problem. To 
derive an effective Hamiltonian describing magnon exci- 
tations, we assume that orbital-spin separation occurs at 
low energies. This is justified when the orbital gap in- 
duced by Ti-O-Ti bond distortions (sec Sec. V) is larger 
than magnon energy. Dynamical coupling between the 
spin and orbital degrees of freedom via fluctuations of 
superexchange bonds and also via on-site spin-orbital in- 
teraction Hso is then considered as a high-energy process, 
leading to an effective spin Hamiltonian. The parameters 
of such a Hamiltonian are obtained by integrating out 
high-energy orbital fluctuations. 

A. Isotropic spin exchange 

We start with estimation of coupling constant J in the 
isotropic spin exchange term, J{Si-Sj). As a flrst step, let 
us consider mean-field approximation, in which the spin 



exchange is given by an expectation value of the orbital 
operator in Eq. (2). Neglecting a small term ''g"^''" {B^^) 

and noticing that (^Ij'') = f-^'Oi one obtains 

^0 = [^\w2 - |(1 -??r2)|£o||r-iJsE. (68) 

The first term (~ 77) is driven by a conventional Hund's 
coupling, while the second one originates from orbital 
singlet correlations in the ground state. For r/ = 0.12, 
these two (classical and quantum) contributions are of 
the same order and give together Jq ^ — 0.214(ri J^e). 

However, the actual value of J measured experimen- 
tally from magnon spectra could in fact be strongly re- 
duced from Jo in Eq. (68) due to a fluctuation effects. In- 
deed, AF and F states are strongly competing in t2g sys- 
tems, and large-scale orbital fluctuations are expected to 
bring about AF spin exchange contribution. We, there- 
fore, have to consider effects of the dynamical spin-orbital 
interaction: 

Hint = Y.^{Si-Sj)6Jij . (69) 

m 

In a ferromagnetic state 5{Si -Sj) ~ — s])(sj 

with sJ being a magnon creation operator. Neglecting 
small -^^rp and ■^-^f^ terms in Eq.(2), exchange inte- 
gral fluctuations are given by SJij ~ ^^i]^ units 
of tiJse)- As the coupling constant in Eq. (69) is not 
small, and because of spins and orbitals may form bound 
states^ in an excited AF states, we will discuss here only 
a qualitative picture. We introduce a correlation func- 
tion D^[^\t) — {TrSJij{0)SJij{T)) describing fluctuations 
of the spin exchange integral. We assume that its spec- 
tral function p(w) = ^D'/-{u) -\- iS) is distributed over 
the characteristic energies larger than low-energy coher- 
ent magnons observed in the experiment (an adiabatic 
approximation which is valid as far as one is concerned 
with low-energy spin excitations). Within this approx- 
imation and neglecting vertex corrections we may eval- 
uate the magnon scattering process on a given bond as 
described in Fig. 9. The result implies a rcnormalization 
of the coupling constant in low energy spin Hamiltonian 
by 6. Jeff = iAj(O) = p{lu)^, which is of AF sign 
as expected. It is the renormalized exchange coupling 
J = Jo + SJef f that determines magnon spectra. We can 
estimate Dij (u) by keeping in SJij ~ 6Aij [Eq. (23)] the 
orbiton pair excitation terms only: 

6jif (pair) = 1 (ajat + bib] -alb]- bid] + H.c). (70) 

We expect that orbiton pair fluctuations are rather inco- 
herent and local, and we parametrize their spectral func- 
tion by a characteristic energy ^pair, obtaining 
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FIG. 9. Scattering of low energy magnons on local fluctusr 
tions of the spin exchange integral, 5Jij . Its Green's function 
(wavy line) is taken at zero frequency (adiabatic approxima- 
tion, sec text). This results in an additional effective spin 
exchange constant (denoted by a filled circle on right hand 
side), 5 Jeff = ^Dii{^) fo"^ energy magnons. 



thus 5 J, 



We notice that Eq. (71) means also 

that {(SJij)^) = |, which can be simply understood as 
follows: on a given exchange bond, e. g., along the c- 
axis, one may have in general nine orbital configurations. 
^From Eq. (11) one observes that these are the orbital 
singlet giving ferromagnetic spin exchange J = — 1, or- 
bital triplet states with J = +1, and five states with 
n\'^^n^^^ = giving zero J. As orbital order is weak, 
and because of the formation of orbital singlets/ triplets 
at given bond necessarily frustrates other neighboring 
bonds, all these configurations will be present giving 
((i5Jy)^) ~ |. It is also natural to expect that flpair ~ 
Worb, with Worb ^ 2 {viJse) sls cstimatcd in the preced- 
ing section. This gives an estimation SJgff ~ | (in units 
of ri Jsb), resulting finally in efii'ective exchange coupling 
for low energy spin excitations as follows: 



J 



-lr]r2-'^il-Vr2)\Eo\ + '^WiJsE- (72) 



These qualitative estimations arc siibstantiatcd in Ap- 
pendix D, in which we calculate magnon energy renor- 
malization within a linear orbital wave theory. 

It is observed from Eq. (72) that J is actually positive 
(antiferromagnetic) for realistic values of ^ 0.12 — 0.13, 
in agreement with the conclusion obtained above from 
energy considerations: The ground state of the model 
is not ferromagnetic in an ideal cubic lattice. At the 
presence of Ti-O-Ti bond angle distortion, J is however 
modified as follows: Eq. (72) obtains a prefactor cos-^ 9, 
and in addition a term — J' ~ — |?7sin^^ [Eq. (58)] has 
to be accounted for. As a result, a classical Hund's rule 
part of J remains unchanged, and the net result 



J 



2 (22 
-^Vr2 + cos^ 0| - - -(1 - Tjra)!^, 
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nJfl (73) 



-0.03ri4°^ 



gives a small ferromagnetic coupling J ~ — u.uj/ 1^'^^ 
for YTiOs with 6 fa 142 deg. Comparing this result (at 
rj = 0.12) with experimental one 

Jexp. 



we obtain the overall energy scale riJ. 



(0) 

SE 



-2.75 meV,i° 

92 meV and 



Jg*^ = 4^0/?/ ~ 59 meV. Bond distortion effect reduces 
the energy scale to riJsE ^ 78 meV and ^ 57 meV in La 
and Y based titanates, respectively. Based on the above 



considerations, we consider J^^ ^ 60 mcV and riJgs 
60 meV as representative energy scales for YTiOs. 

A main message of the above considerations is that 
the spin-exchange constant as seen by a coherent low- 
energy magnon excitations in YTiOs represents in fact 
only a small fraction of the real strength of dynamical 
spin couplings. Because the sign of f2g-spin exchange is 
not unique, and because the orbital order is weak, large 
fluctuations of the spin couplings are present in titanates. 



B. Anisotropic SE interaction 

Next, we consider effects of a relativistic spin-orbit cou- 
pling 



(74) 



This interaction introduces anisotropy in the effective 
spin Hamiltonian, which obtains (besides the rotationally 
invariant Heisenberg part) an additional, so-called anti- 
symmetric Dzyaloshinskii-Moriya (DM) and symmetric 
anisotropy interactions.^^ The anisotropic interactions 
select orientation of the magnetization in the crystal, and 
lead also to magnon gap(s). The structure of anisotropic 
terms is essentially determined by orbital state via ex- 
pectation values and dynamics of the angular momentum 
operator in Eq.(74). Thus, we would like to obtain spin- 
orbit coupling induced corrections to the spin Hamilto- 
nian, and discuss their consequences on magnon spectra, 
thereby testing the proposed orbital state. 

As usual, anisotropic interactions are obtained by per- 
turbation theory involving both isotropic Hse and Hso- 
We mostly discuss the quadrupole ordered orbital state 
with condensed c orbitals (see Sec. IIIB). In derivation 
of the anisotropy Hamiltonian, we need to keep in the 
superexchange operators J^j ^ and kIJ'' in Eq.(l) such 

terms that (i) operate in the ab excited states and/or (ii) 
connect a ground state with excited states of orbitals. 

For instance. 



7(c) 



JSE 



ia + rija + aid] + didj - ajftt - dibj 



+{a^b)]. 



(75) 



[Here, the terms proportional to small numbers (ri — r2) 
and (r2 — r^) are neglected]. 

We consider nearest- neighboring sites i and j. The 
local excitation energy to create an orbiton a or is de- 
noted as Aloe- It is reasonable to associate Ajoc with the 
"center of gravity" of the orbiton band that covers the 
energy window from ^ ri Js^; to ^ 2 riJsE as obtained 
in Sec. V. Thus, we will consider A(oc ~ 1.5ri Jsb in our 
estimations, when we compare later on the results with 
experiment. 

To obtain spin anisotropy interactions, it is convenient 
to work again in a rotated basis, applying transformation 
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(19) also for the spins. The scalar product of NN spins 
in HsE is then expressed as Si ■ Sj = SiT^'^^Sj where 



rp{a) 




(76) 



and 



, / -(It 2s) -1 l±2s 

T'-^'^^ = - I -1 -(l±2s) 1=f2s I (77) 

^ \ l±2s 1=f2s -1 

for the state I (a). [For the state I (b), matrices T for 
bonds a and b are equal and given by Eq. (76)]. 

Third-order perturbation with respect to Hgo and Hse 
gives a symmetric part of the spin anisotropy Hamilto- 
nian Hani in a local coordinate: 



with 



H. 



M('=) = - 
2 



(7) 



1 c-3s 

2 c + 3s 



(78) 



(79) 



c — 3s c + 3s 



-1 



Anisotropy constant A in Eq. (78) is given by the follow- 
ing expression: 



4 rv+rs/ A 



(80) 



Transformation of the spin operator from a local to global 

coordinates is expressed as Si = RiRSi. The matrices Ri 
{i = 1,2,3 and 4) transforming spin coordinates at four 
sublattices are as follows: 






I -I 0\ 


)■-( 


0-10 










j ,^4 = 






U -1/ 



(81) 



Using this transformation, one obtains a symmetric 
anisotropy Hamiltonian defined in the global coordinate: 



H^^ = -ASiMijSj, 



(82) 



where depends on the NN bond. For 1-3 (2-4) bonds 
along the c axis, Mi3(24) is given by 



/3/2 ±3 
Mi3(24) = 7 7/2 
* \ ±3 -3/2 



(83) 



The interaction matrices for NN spins on a, b bonds have 
similar structure: 



-yiyj 



.1a 



ytyj-^i^j 



-> a 



FIG. 10. Schematic representation of the structure of sym- 
metric spin anisotropy interactions of e^-symmetry. Inter- 
actions along different bonds are denoted by OiOj, which 
should read as SiaSja times an overall interaction constant 
given below each figure. For instance, a-bond interactions are 
— lASizSj:, and ^A{SixSjx — SiySjy) (note eg(3z^ - r^)- and 
eg{x^ — 3/^)-type symmetry), where the constant A is defined 
by Eq.(80). Overall cubic symmetry of the interactions is 
evident for both contributions. 



M- 



12(34) 



M- 



14(23) 



-3/2 ±3 
±3 3/2 
7/2 

7/2 
-3/2 ±3 
±3 3/2 



(84) 



Symmetric anisotropy interactions can be classified ac- 
cording to cubic invariants: there are terms of Cg and 
t2g symmetries, generated by diagonal and nondiagonal 
elements of the matrices Alij , correspondingly. For con- 
venience, we show the interactions for the state I (a) in 
Fig. 10 and Fig. 11 (a), in which bond dependence, direc- 
tion of spins and the scale of individual anisotropy terms 
are shown. The Cg symmetry anisotropy in Fig. 10 (a) 
has been discussed in Ref. 13 under the name of "cubic" 
anisotropy in the context of magnon gap in LaTiOs. A 
remarkable feature of this interaction is its intrinsic frus- 
trations: namely, treated classically, it acquires a rota- 
tional symmetry in spin sector, resulting in an infinite 
degeneracy of classical states. An accidental pseudo- 
Goldstone mode, which appears in classical limit, can 
acquire finite gap by quantum fluctuations only. In the 
present orbital ordered states, the symmetry is lowered 
compared with the orbital liquid state in LaTiOs. Thus, 
additional terms are generated as shown in Figs. 10 (b) 
and 11 (a). It is noticed that these terms have a similar 
frustrated nature: summed over all the bonds, they can- 
cel each other exactly. Thus, only a small gap is expected 
from these interactions. 

Physically, the structure of anisotropy interactions is 
determined by local correlations of the angular momen- 
tum, and can therefore be traced back to noncoUinear 
arrangements of these correlations shown in Fig. 2 (a). 
For instance, the leading, "cubic" term [see Fig. 10 (a)] 
reflects that ly components are correlated fcrromagncti- 
cally along the c axis, while Izilx) components are parallel 
along a (6) axes. 
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FIG. 11. Left: Schematic representation of the structure of 
symmetric spin anisotropy interactions of t2g-symmetry. The 
notation {af3)ij stands for SiaSjfj + SifjSjoi multipUed by in- 
teraction constant "iA/A. Black arrows represent the direction 
of spins ftivoral)lc for this interaction when the direction of the 
uniform moment is taken along [001]. Right: Antisymmetric 
DM spin anisotropy interactions [Eq.(88)]. Gray arrow de- 
noted by dij shows the orientation of DM vectors on different 
bonds. A preferred spin pattern for this interaction is shown 
by black arrows. 



In the state I (b) , the leading anisotropy interaction in 
the a6-plane is given by —S^S^ and that along c-direction 
is obtained to be —SfSj [as can easily be seen also from 
Fig. 2 (b)]. The number of bonds with anisotropic inter- 
action —SfSj is two times larger than that with —SfSj. 
This breaks the rotational symmetry in a spin sector even 
in the classical limit, and generates a large magnon gap 

Antisymmetric DM interaction appears in a second- 
order perturbation theory as a combined effect of Hse 
and Hgo- The remarkable feature of the orbital state I is 
that leading terms in SE interactions, which arc propor- 
tional to (ri -t-r2) in Eq.(2), do not contribute to the DM 
interaction. That is because of the classical expectation 
value of Aij in Eqs. (4) and (8) vanishes in our orbital 
states. Rather much smaller Hund's coupling terms pro- 
portional to the small number (n — r2) only give rise to 
DM interaction. This feature contrasts with that in the 
orbital state reported in Refs. 15, 23-26, in which a large 
DM interaction is present (see Appendix E). 

After somewhat tedious but straightforward calcula- 
tions one obtains the following interaction between NN 
sites: 



DSiN^^^Sj. 



Here, the matrix iV^*^^ reads as 
1 



3 



2c - 2s 1 c+s 

1 2c+2s c - s 
c + s c — s —2 



(85) 



(86) 



and the interaction constant D is obtained as follows: 
ri - r2 A 



Transforming the local spin axes to the global ones, one 
arrives at the following DM interaction: 



H 



(ij) 
DM 



Ddii 



(88) 



with dij = a-. Here, dj is the unit vector parallel to one 
of local axes (.x^, i/; and Zi) which is perpendicular to 
the i-j bond direction and antiparallcl to its counterpart 
at site j [see Fig. 2 (a)]. For example, rfia = (1,0,0), 
di2 = (0,1,0), fii4 = (0,0,1), etc. For convenience, we 
show the DM interaction for state I (a) in Fig. 11 (b). 



C. Spin waves: comparison with experiment 

Now, wc discuss the anisotropic spin interactions in 
the context of the experimental observations of magnon 
dispersion in YTiOs.^" 

Cubic symm,etry of the spin wave dispersion: This 
puzzling observation is naturally explained by the present 
theory, as ferromagnetic couplings in all the states I and 
II are perfectly isotropic, = j(^) = The reason 
is high symmetry of the orbital ordering patterns, as can 
be visualized from Fig. 4. It is stressed that this result 
is robust, "no fine-tuning" property of the model. (The 
isotropy is expected to be relaxed somewhat by lattice 
distortions. However, effects of two types of distortions, 
that is, Ti-O-Ti bond angle distortion and elongation of 
TiOe octahedron, on the anisotropy of spin couplings are 
opposite and almost cancel each other.^^) 

Isotropy in spin space, magnon gap: 
Magnons in YTiOs are almost gapless^" (upper limit for 
the gap is 0.3 meV); this is a serious test for possible 
orbital orderings. We show now that high symmetry 
of orbital orderings in the present theory resolves this 
problem as well. The crucial point is the frustrating na- 
ture of the anisotropic spin interactions obtained above: 
Even though anisotropic couplings on an individual bond 
are substantial, there is cancellation of classical contri- 
butions stemming from different bonds, and one obtains 
only small gaps of quantum origin. We illustrate this by 
considering first the leading "cubic" term. 

Quantum magnon gap by "cubic" term: Consider the 
effective spin Hamiltonian in the state I (a) given as fol- 
lows: 



H>i = Hi! 



Ha 



(7) 0(7) 



Y,[JSi-Sj+ASr'S^ 

{ij) 



(89) 



D = JsE- 



6 A,. 



(87) 



Here, Hiso represents the isotropic SE Hamiltonian. Co- 
efficient of the "cubic" term is A = 7/8A [see Fig. 10 (a)], 
and axes in the spin space are changed as S^, 
and S*^ such that 7 corresponds to the 

direction of the i-j pair. Thus, we have —SfSj for c- 
bonds, —SfSj for a-bonds and —SfSj for 6-bonds. 
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Because the effective spin Hamiltonian Hs has a dis- 
crete (cubic) symmetry, magnon excitation is expected 
to have a gap. However, due to the rotational sym- 
metry of Eq. (89) in the hmit of classical spins, linear 
spin- wave theory cannot provide finite gap. This prob- 
lem is resolved by the order-from-disorder mechanism, 
which selects a particular classical state which provides 
the largest zero-point energy when fluctuations are in- 
cluded. This opens also a magnon gap.^'^'^*^ Thus, we 
calculate spin-wave contribution to the ground-state en- 
ergy as a function of the angle 6 between the c axis and 
the uniform moment. First, we rotate the spin quantiza- 
tion axes around the b axis: 

Sf = cgS^ + seSl, 

S! = -S8Sf + ceS^, (90) 

where ce = cos 6 and sg = sinO. Second, by us- 
ing Holstein-Primakoff approximation, we obtain the 
magnon dispersion which shows explicit 6 dependence as 

LOj^ie) = zSJ^Xj:[Yj:^a{c,-c^)s% (91) 

where = (l-7i) + a(l-Cj,), Y^, = (1 - 71) + a(l - c^), 
and a = A/ 2, J represents the ratio between anisotropic 
and isotropic interactions. Finally, by calculating the 
zero point magnon energy, we obtain the ground-state 
energy Eq per site. In the limit of 6' ^ 1, 

E^{e) = —zJSiS + 1)(1 + a) + l^ a;g(e) 

k 

= -const + K^ffS'^e'^. (92) 

Here, Kgff represents an effective spring constant. In 
the case of a 4C 1 which we are interested in, the spring 
constant is given by -K'e// = A'^R/zSJ with 

2 

i? = 3^ « 0.28 . (93) 

k 

The potential Eo{9) in Eq.(92) can be associated with an 
effective uniaxial spin anisotropy Hamiltonian 

H^f} = -Keff E St SI (94) 

generated in the symmetry-broken phase (with spins ori- 
ented along [001]). Therefore, one finds a magnon gap 
^mag = '^SK^ff = 2A^R/zJ. Note that ls.mag is inde- 
pendent of S and is proportional to A?. We confirmed 
that the same magnitude of the magnon gap is derived 
from Hs using single mode approximation^^ (see Ap- 
pendix C). Thus, "cubic" anisotropy gives the magnon 

gap I^mag = 2{lAfR/zJ. 

Now, we estimate the coupling constant A. ^From 
Eq. (80) with nJsE ~ 60 meV and r? ~ 0.12, one 



finds A ~ 23(A/A;oc)^- As discussed in the precedings 
section, we consider Ajoc ^ l-S^iJgg ^ 90 mcV. Us- 
ing the atomic value Aat = 19 meV,^° one obtains then 
A ~ 1 meV consistent with the experimental value of 
A obtained in Ref. 10 from magnon spectra. With this 
value of A, the magnon gap from the "cubic" term is only 
K, 0.03 meV in the state I (a). 

On the other hand, a large classical magnon gap 
Amag — jV^SA w 1.2 mcV is obtained in the state 
I (b), as we already discussed in the preceding section. 
This allows us to exclude the orbital ordering structure 
I (b) from possible candidates for YTiOa, although we 
do not know precisely which kind of lattice distortions 
favor state I (a) over configuration I (b). 

Contributions from the other terms: Anisotropy of 
eg{x'^ — y'^) symmetry [see Fig. 10 (b)] can be analyzed 
similarly; we find that it also supports an easy magne- 
tization axis along one of the cubic axes, say [001]. Its 
contribution to the magnon gap, A^^g = 6{^A)^R/zJ ~ 
0.014 mcV is smaller than that of the "cubic" term, as 
expected. 

Once the [001] direction is chosen as the direction of 

uniform moment, spin anisotropy terms of t2g symme- 
try and DM interactions give rise to spin cantings as 
shown in Figs. 11 (a,b). The canting angle 6 is given 
by 61 w (3/2\/2)(^/4J) ~ 0.1 rad for t2g symmetric 
anisotropy interaction, and 9 ~ ^/2{D/AJ) ~ 0.07 rad 
for the antisymmetric DM interaction. [DM interac- 
tion constant D « 0.57 meV is estimated from Eq.(87)]. 
These values are within the experimental canting an- 
gles ~ 0.17 rad.^° Finally, contributions of t2g symmet- 
ric anisotropy and DM interactions to the magnon gap 
Amag are estimated as « AJS6'^{t2g) ^ 0.05 meV and 
« AJS9'^{DM) ^ 0.03 meV, correspondingly; these num- 
bers are rather small again. It should be noticed, that a 
more quantitative analysis of the problem, in particular 
the precise structure of the spin canting pattern requires 
consideration of all the anisotropy terms on equal footing, 
which will be presented elsewhere. 

To summarize this section, we have obtained spin 
anisotropy interactions induced by the spin-orbit cou- 
pling, and considered their effects on spin-wave spectra. 
Because of high symmetry of the orbital ordering, par- 
ticularly in the quadrupole ordered state I (a), magnon 
dispersion is found to have cubic symmetry, magnon gap 
is small, and there are small cantings of spins away from 
the c axis. All these observations are consistent with 
experiment. 

VII. ORBITAL CONTRIBUTION TO THE 
RESONANT X-RAY SCATTERING 

We turn to the discussion of further experiments which 
may help to verify the proposed orbital state in YTiOs. 
First, we consider the resonant x-ray scattering, which 
has proven to be a useful method in the study of or- 
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bital order symmetry.''^ ■"'^ The following section will be 
devoted to possible ways of detecting orbital excitations. 

We focus on the orbital state I (a) [which is the most 
plausible candidate as discussed in previous sections]. 
While the exchange bonds in this state are the same 
(important for the isotropy of spin waves), a local sym- 
metry is lower than a cubic one [see Fig. (4)]. Thus, or- 
bital order may induce spatial modulations of the level 
structure of an excited photoelectron in p-states via the 
so-called Coulomb mechanism. This may lead to addi- 
tional weak reflections at orbital ordering vectors. Pre- 
dictions of our theory for such an experiment are as fol- 
lows. Orbital order shown in Fig. 4(1) is identified as a 
three-component quadrupole ordering of t2g-symmetry, 
Ta{q) with a = x,y,z. Each component has its own 
propagation vector: 

Tx = {lylz + ^zly) ^ = 3*^^ ' 

Tz = {Uv + lylz)R = \Qe'^'-^, (95) 

where qi = (7r,0, tt), q2 = (7r,7r,0), and = (0,7r,7r). 
As order parameter Q is strongly suppressed by quan- 
tum fluctuations {Q ~ 0.19, see Sec. Ill B), we obtain 
that each component has only a small amplitude, giv- 
ing \Ta\^ ~ 0.016. This imphes that the corresponding 
anomalous Bragg intensity is at least 60 times weaker 
compared with the classical orbital orderings. It might 
therefore be very difScult to single out this contribution. 
However, new azimuthal (</?) and scattering [Og] angu- 
lar dependencies of an additional intensity, which should 
show up below Torb, may help to identify order symmetry. 
For (7r,7r,0) [(100) in orthorhombic notations] scattering 
these dependences are obtained as follows 

laa' (tt, tt, 0) DC sin^ 2(/? (96) 

for <j-cf' polarization (see for notations Ref. 31), and 

-^cttt' (ti"! TT, 0) cx (cos2(^sin^s -|- sint/jcos^s)^ (97) 

for (T-tt' polarization. (Azimuthal angle <^ = corre- 
sponds to the configuration in which the diffraction plane 
is parallel to the c axis.) Scattering intensities at qi and 
5*3 (which are contributed by and components, re- 
spectively) can be obtained from symmetry considera- 
tions. 

It is still a controversial issue whether the resonant 
x-ray scattering, observed in YTiOs,^^ is related to or- 
bital order or lattice distortions (see Refs. 33, 26, 17). 
Either way, we expect that the orbital order contribu- 
tion, if present, must be temperature dependent reflect- 
ing orbital order /disorder transition, as in the case of 
manganites with strong orbital order. Therefore, a care- 
ful analysis of the T-dependence of reflections at orbital 
ordering vectors is desirable. 



Above discussion brings us to the problem of the or- 
bital ordering temperature in YTiOs. Thus far, there are 
no reports on the orbital ordering temperature in YTiOs 
(weak structural change at spin ordering temperature^^ is 
only an indirect indication). In our SE- model picture, we 
expect that this transition should occur at low tempera- 
tures only, and we suspect in fact that orbitals in YTiOs 
do order at ferromagnetic spin transition Tc- This is be- 
cause of strong spin/orbital coupling, and also because of 
the frustrated nature of orbital-only model itself, which, 
as we have shown, may develop long-range order on a cu- 
bic lattice at zero temperature only. Lattice distortions 
that open a finite orbital gap allow finite-temperature 
transition, but this cannot occur much above the fer- 
romagnetic transition, because of strong disorder intro- 
duced by spin fluctuations in the paramagnetic phase. In 
other words, orbital order and isotropic spin ferromag- 
netism are intimately connected, supporting each other. 
Physically, this implies that short-range ferromagnetic 
correlations are of vital importance for orbital ordering, 
and vice versa. 

A quantitative description of the finite-temperature 
behavior of a realistic spin-orbital model is complicated. 
We may give only very rough estimation for the orbital 
ordering temperature based on the mean-field picture. As 
we are going to ignore fluctuations completely, this esti- 
mation should be regarded as an upper limit, which we 
would like to know. To this end, we consider a spin para- 
magnetic phase and set {Si ■ Sj) = in Eq. (1), neglect in 

the orbital interactions ^l^j' in Eq. (23) all the terms ex- 
cept those which contain an emerging quadrupole order 
parameter. This leads to 

Horb ~ TT ^ ^ QizQjz 
oZ 

{^3) 

=^-l{Q)Qiz- (98) 

From Eq. (98) we obtain Torb = g (in units of riJss), 
which, including the Ti-O-Ti bond angle {6 ~ 142 deg) 
correction for YTiOs, reads as 

Torb = I cos2 0(ri4"]) ~ 0.1 (ri4°i) . (99) 

On the other hand, spin ordering (mean-fleld) temper- 
ature is Tc = ^\J\, with J ~ — 0.03(ri Jg^) given in 
Eq. (73). Both Torb and Tc should, of course, be reduced 
by fluctuations (indeed, with J^^p = —2.75 meV,^° one 
obtains mean-field — 48 K instead of observed 27 K), 
so it makes more sense to consider their ratio, which is 

Torb/Tc ^ 2.2 . (100) 

For Tc ~ 25 — 30 K (which is sample dependent), this 

gives an upper estimation Torb ^ 55 — 70 K. We would 
like to think that local orbital order, accompanied by 
short-range ferromagnetic correlations, starts to develop 
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at these temperatures. In fact, the presence of such a 
correlations in YTiOs below ~ 50 — 60 K has been re- 
ported from several experiments: (i) sharp drop in NMR 
relaxation rate, which has been speculated in terms of 
orbital ordering,^'* (ii) spin-resonance line shape changes 
from typical paramagnetic spectra to the ferromagnetic 
one,^^ (iii) weak quasielastic magnetic scattering is ob- 
served above Tc^° 



VIII. ANGULAR MOMENTUM FLUCTUATIONS: 
DYNAMICAL MAGNETIC SUSCEPTIBILITY 

In this section, we would like to calculate orbital contri- 
bution to the inelastic neutron-scattering intensity. The 
point is that t2g orbitals are magnetically active, as their 
angular momentum may directly couple to the neutrons. 
Of course, there is a contribution also in nonmagnetic 
channels: The orbital quadrupole moment is coupled to 
the phonons, and hence single or double orbiton (depend- 
ing on the structure of this coupling) may be excited by 
neutrons indirectly via lattice vibrations. We focus on 
the magnetic scattering, and calculate orbital angular 
momentum dynamical susceptibility. If the t2g orbital 
level is split up by strong lattice distortions, one would 
expect just a local, crystal-field transitions. In SE-driven 
orbital picture, advocated in this paper, angular momen- 
tum fluctuations are however of the collective nature. 
Thus, we expect momentum selected (though strongly 
damped) transitions, forming broad bands. 



A. Quadrupole order 

Consider first local angular momentum susceptibility 
Xzoc(w) in the quadrupole ordered state. It is defined as 



(101) 



and its imaginary part describes the spectral shape of the 
momentum-integrated inelastic neutron-scattering cross 
section. 

In a linear orbital wave approximation, the imaginary 
part of xioci'^) at a; > is given by 



In order to account for a finite gap induced in the orbital 
sector by lattice distortions (see Sec.V), we use hereafter 
Eq. (66) for the orbital excitation spectrum. The numer- 
ical result for x'l'od^) is presented in Fig. 12. The sharp 
structure about Worb (taken as an energy scale in the 
figure) is related to the orbiton band-edge effects, which 
should go away when damping effects are properly taken 
into account. 
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FIG. 12. Solid line: Imaginary part of the local angular 
momentum susceptibility Xiod'^) [Eq.(102)] in quadrupole or- 
dered state of orbitals. The energy uj is given in units of 
Worb defined by Eq.(67). A finite gap for the orbital waves 
stemming from symmetry breaking interactions is taken into 
account according to Eq. (66) with parameters / = 0.086 
and e = 0.18 (this gives Worb — 2riJsE)- The sharp peak 
structure is expected to be smoothed by damping effects (not 
accounted in the present study), as indicated by the broken 
line. 



We turn now to the momentum dependence of dynam- 
ical susceptibihty, x((f, uj) = (j,^ ■ . This quantity 

determines a dynamical structure factor, which at T = 
is given by 



S{q,uj) = -Imx{q,uj). 

TT 



(103) 



A noncollinear, four sublattice orbital order leads to the 
following structure for x(g,a;): 



x(9,w) = -[xiiq + qi) + xi{q + Q2) + xiiq + Q3)]^ 

+ ^[X2{q + q2)-X2{q + qi)]^ (104) 

+\ [2x3(9 + 93) - xaiq + qi) - xaiq + 9*2)] ^• 



Here, the orbital ordering vectors 



qi = (7r,0,7r), 
0*2 = (tt, 7r,0), 
93 = (0,7r,7r) 



(105) 



(102) for the state (a) and 



qi = (7r,7r,7r), 

02 = (7r,7r,0), 

03 = (0,0,7r) 

for the state (b). The susceptibilities 

Xi{q,u^) = \{W%^+lpl^)^, 



(106) 
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(107) 



arc defined in a rotated basis given by the transforma- 
tion in Eq. (19). We calculate these susceptibilities in a 
linear orbital-wave approximation. The imaginary part 
of Xa{q, cli) at cj > is obtained as 

+b J^^'~'' 5{u-U2^), (108) 
where Ai = Bi = 1/2, A2 = -B2 = 72/2K and A3 = 

-B3 = 73/2K. 

Numerical results for S{q,uj) in the state I (a) are 
shown in Fig. 13. An intensive hot spot at momentum 
q = (tt, tt, tt) at energies about orbiton gap, and flat dis- 
persions at (tt, TT, g^) direction are noticed. The rather 
complicated multiband structure has its origin in non- 
collinear nature of the underlying orbital ordering, char- 
acterized by a several ordering vectors. [In the state I (b), 
which has different ordering vectors qi, S{q, lo) shows dif- 
ferent energy- momentum structure (not shown)]. The 
energy scale for orbital fluctuations (~ Worb ^ 2ri Jg^) 
is much larger than magnon energies. This is because 
of strong cancellation of ferromagnetic and AF contri- 
butions to the spin-exchange integral J (see Sec. VIA), 
resulting in rather small magnon bandwidth (which is 
only a fraction of tiJse)- Therefore, magnon excitations 
are expected to be well defined, since they are located 
within the orbital gap. As for the high-energy orbital 
excitations, we expect strong damping effects stemming 
from nonlinear couplings between orbital waves them- 
selves, and also from the dynamical coupling between 
spin and orbital fluctuations. These effects should in 
fact relax momentum resolution and smooth away sharp 
structures obtained in Fig. 13 by using undamped orbital 
waves. 



B. Orbital magnetic state 

For completeness, we also give equations for the mag- 
netic response of the magnetically ordered state of or- 
bitals. In contrast to the quadrupole ordering, this state 
gives rise to static Bragg peaks of orbital origin. These 
peaks are located at orbital ordering vectors qi [Eqs.(105) 
and (106)]: 

(^Iq-^-q) = ^("^') [^(^^ ^1) + S{q- q2) + S{q~ q^)], 

(109) 

and their intensity is determined by the orbital magnetic 
order parameter mi — 0.19 [Eq.(38)]. 
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FIG. 13. (Color online). Intensity of the orbital contribu- 
tion to the magnetic structure factor S{q, u>) [Eq.(103)] in the 
quadrupole ordered state 1(a). Energy ui is given in units of 
Worb defined by Eq.(67). A finite orbital excitation gap due 
to symmetry breaking terms (see for details Sec.V) is taken 
into account according to Eq. (66) with parameters / — 0.086 
and e = 0.18. 



The dynamical susceptibility x((f, w) is given by the 
same form as Eq. (104) with Xi.2,3('?, i-^) having the same 
definitions given by Eqs. (107). In the orbital magnetic 
state, these susceptibilities are obtained as follows: 



x'l{q,^)^j 



— 6{UJ - LUl^) + —6{UJ - UJ2q) 



(110) 



and 



xLC^,^) 



7''72,3 

An 



d{UJ - LOl^) d[lAj - UJ2^) 

^Iq ^1q 



IX. SUMMARY AND DISCUSSION 

In this paper we investigated a spin-orbital superex- 
change Hamiltonian in a Mott insulator with elec- 
tron configuration, focusing mainly on the orbital order 
and dynamics in the spin ferromagnetic state. An im- 
portant feature of the Hamiltonian in the spin polarized 
state is the large frustration of orbital states, thus the 
ground state is governed by the interplay between or- 
bital frustration and quantum fluctuations. On the clas- 
sical level, there is a local Z2 symmetry which leads to 
an infinite degeneracy of classical configurations. Long- 
range orbital order does occur in the model by a quantum 
order-from-disorder mechanism, which selects a particu- 
lar ordering patterns. Orbital orderings are quite unusual 
having highly noncoUinear four sublattice structure, and 
provide the same spin couplings in all cubic directions. 

Besides classical local Z2 symmetry which is removed 
by quantum dynamics, there are exact conservation laws 
in the orbital Hamiltonian. They are related to the con- 
servation of orbital quantum numbers in the SE process, 
and lead to a multitude of degenerate quantum ground 
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states which can smoothly be connected to each other 
by phase rotations of the complex orbital order param- 
eter. Such continuous rotations generate orbital Gold- 
stone modes, which have 2D dispersion because of planar 
geometry of t2g orbitals. As a result, static orbital order 
sets in at zero temperature only. Degenerate quantum 
ground states are physically different: depending on the 
phase of the orbital condensate, they describe quadrupole 
or magnetic orderings or their coherent mixture. Extrin- 
sic perturbations, e.g. lattice distortions or spin-orbit in- 
teractions may remove the degeneracy and fix the phase 
of the condensate. Reflecting the large quantum fluctua- 
tions, the orbital order parameter is unusually small. 

We found that the orbitally ordered ferromagnetic 
state is slightly higher in energy than the spin-AF or- 
bital liquid state. This is because the latter state gains 
an additional quantum energy from coupled spin-orbital 
fluctuations. To explain ferromagnetism of YTiOs, we 
emphasized the role played by Ti-O-Ti bond angle dis- 
tortion. This distortion favors ferrostate by generating 
an unfrustrated ferromagnetic SE interaction via virtual 
hopping of electrons between NN t2g and Cg orbitals. 
Even more importantly, the bond distortion eliminates 
orbital soft modes and opens a large orbital gap, such 
that orbital order becomes stable at finite temperature. 
This distortion stabilizes the quadrupole ordered state. 

The strong competition between AF to F states in 
the present model has direct relevance to nearly continu- 
ous transition between these states observed in titanates. 
In these compounds, A-site substitution from La to Y 
increases the Ti-O-Ti bond distortion, hence changing 
gradually a delicate balance between AF and F states. 
Because of the orbital fluctuations, spin-exchange inte- 
gral on every link experiences strong fluctuations, both in 
amplitude and in sign, and the system may develop either 
an AF or F state depending on local orbital correlations. 
In this picture of "fluctuating exchange bonds," the mag- 
netic transition temperatures, Tpf and Tq, represent only 
a time-averaged static component of the spin couplings. 
Its value is only a fraction of full superexchange energy 
scale, and can gradually be tuned by external forces such 
as lattice distortion, pressure, etc. We think that weak 
orbital order may continuously evolve in titanates when 
the bond angle decreases below a certain critical value, 
and propose the phase diagram shown in Fig. 14. The 
sign of the time-averaged spin coupling depends on a 
local correlation of orbitals. To the right of the criti- 
cal point orbital correlations are more antiferromagnetic, 
supported by noncollinear orbital orderings. To the left, 
the genuine ground state of the t2g superexchange, an 
orbital disordered state supporting spin AF is stabilized. 
In the proximity fluctuating part of the overall 

superexchange interaction dominates, and separation of 
the spin and orbital degrees of freedom might no longer 
be possible. This scenario can be tested experimentally 
by investigating the spin and orbital transition temper- 
atures under high pressure and magnetic field. On the 
theoretical side, a quantitative description of the transi- 
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FIG. 14. Proposed picture for the evolution of magnetic 
and orbital states in perovskite titanates. Tc and Tjv arc 
Curie and Neel spin ordering temperatures, respectively. Or- 
bital ordering below Torb — Tc is expected in ferromag- 
netic region. Below T*, short range ferro-type spin and non- 
collinear, dominantly "AF"-type orbital correlations grow up. 

tion between ferromagnetic and AF states, controlled by 
orbital order-disorder transition, remains an interesting 
and challenging problem. 

Further, in order to discuss the recent spin-wave data 
in YTiOa, we derived the low-energy spin Hamiltonian by 
including a relativistic spin-orbit coupling and lattice dis- 
tortions that induce the orbiton gap. Using this Hamil- 
tonian, we calculated the magnon gap and found that it 
is very small because of high symmetry of the underlying 
orbital order. Both the real and the spin space isotropy of 
the spin-wave spectra observed experimentally find nat- 
ural explanations within the proposed theory. 

We also calculated the orbital contribution to resonant 
x-ray intensity and to dynamical magnetic structure fac- 
tor. Predictions made should be helpful in further exper- 
imental study of titanates. 

ACKNOWLEDGMENTS 

We would like to thank B. Keimer and S. Maekawa 
for stimulating discussions. Discussions with C. Ulrich, 
S. Ishihara, T. Kiyama, M. Itoh and J. Akimitsu are also 
acknowledged. One of us (G. Kh.) would like to thank 
International Frontier Center for Advanced Materials at 
the Institute for Materials Research, Tohoku University, 
where a part of this work was carried out, for its kind 
hospitality. S. O. acknowledges the hospitality of Max- 
Planck-Institut fiir Festkorperforschung in Stuttgart dur- 
ing his several visits. 



20 



APPENDIX A: BOND DEPENDENCE OF 
QUADRUPOLE MOMENT OPERATORS 



APPENDIX C: MAGNON GAP BY 
SINGLE-MODE APPROXIMATION 



In this appendix, wc present the cxphcit expressions 
for Q^^\ T^K By replacing /„ in Eqs. (24) with one 
obtains 

(7) _ / -cQx T sTz, for 7 = a{h) 



Qx, for 7 = c 



^(7) ^ / -cT;±sQx, for7 = a(6) 

\ T„ for7 = c 

^(7) ^ f -cTx T sTy, for 7 = a{b) 

\ Tx, for7 = c 

7^(7) = / -cTy ± sT^, for 7 = a{b) 

1^ Ty. for 7 = c ^ ' 

Til^ is given by T^,^ = T^^^ ± T^^\ 



APPENDIX B: COEFFICIENTS OF 
BOGOLIUBOV TRANSFORMATION 

We present the coefficients of Bogoliubov transforma- 
tion which diagonaUze How in Eq. (28). This transfor- 
mation reads as follows: 

= wcosh^itt^j -\- V cosh. 920L^j^ — usinh^iQ^ - 



-V sinh 6201 



t 

2 -fe 



6g = — w coshf^iQ'j^j: + ucosh.920L^^ + tisinh^iaj^ - 
-usinh6i2a^_g. (Bl) 



The inverse transformation is given as follows: 



^ife ~ wcosh^iag — t;cosh^i6jg -|- usinh^ia^ - 



-t; sinh ^16^ -. , 

— k 



'^2k ~ vcosh92aj: + ucosh62bj: -|- i'sinh02a^ j 
-|-wsinh02^^ - • 

—k 



(B2) 



Here, u and v are 



1 L ^12 

u = H 

V2\ K 



(B3) 



We examine here the magnon excitation gap, generated 
by "cubic" anisotropy interaction given in Sec. VI B, by 
using a different approach: namely, we apply single-mode 
approximation.^^ In this approximation, spin excitation 
energy at g ^ is given as follows: 



(Cl) 



where H, is given in Eq. (89) with S^"'^ = fi-f , S'f ^ = Sf , 
and S^'^'^ = Sf. The double-commutator correlation func- 
tion in the numerator of Eq. (CI) is equal to X^^^ ^-J^ 
where 



rib) 



-AiStS'j 



Sj Sj 



T^f = -a( - 2SfS' + StS: 



(C2) 



By using linear spin-wave theory, one obtains 



(ij) k 

= 2S{l + alao}. 



1 



«/g--^0 



(C3) 



Expectation values (aia^) and (ata^ -) are calculated at 

T = 0. Up to linear in A terms, one obtains (aia^) = 
and 



A 



2zJ 1 - 71 
Consequently, we find the magnon gap 

2A'^R 



zJ 



(C4) 



(C5) 



where R is given by Eq. (93). This is exactly the result 
obtained in Sec. VIC. 



and 



1 / 72 



(B4) 



respectively, and ^1,2 is given by tanh2^i^2 = 71 ± k. 



APPENDIX D: EFFECT OF ORBITAL 
EXCITATIONS ON THE MAGNON DISPERSION 

In Sec. VIA, we discussed the renormalization of 
nearest-neighbor isotropic spin coupling J by orbital fluc- 
tuations. Here, we investigate this effect in more detail, 
by considering effects of the dynamical spin/orbital cou- 
pling on magnon spectra. In terms of magnon and 
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orbiton a^, operators, the dynamical spin/orbital cou- 
pling in Eq.(69) is expressed as (in units of riJss) 

Hint = -^^^I'Sp Q^o(ai,a-^+bi,b-^ 



pqq> 



(ri + r2)44+(ri-r2)6p;- 

H.c. (Dl) 



where j7 ^ p — q — (f . Factor Q ~ 0.19 stems from 
Hartree decoupling of the QizQjz term in Eq. (23). Fq 
and Fi(^i.2,3) are given as follows: 



To — 1 + 7l q+q' - I'l 7l q+q'-p 

= 7i (f + 7* (f ^ I'i. q-p ^ I'i q'-p- 



(D2) 



Second-order perturbation with respect to Hint gives 
the renormalization of magnon excitation energy as up = 

u)^^ — 6ujp. Here, magnon softening 5u)j; is given by 



<i<f 



M(l + A) 



+ 



+ (^l? + — ^p 

M(l + A) 



Wig- + U)2q' + ^j? — Wp- 



(D3) 



where 0Ji^2q = {iJ^ — 1±Y^'^ with 7± = 71 ± k correspond 
to two orbiton branches. The orbiton chemical potential 
/U controls the orbital gap. Magnon energy u)p that enters 
in this equation is considered to have a NN Heisenberg 
form 3| J|(l — 7ig-). The matrix elements in the numerator 
have a following structure: 



M 



{Tr+rT,)\l+x)+Tl{l + y) + Tl{l-y) 

+2(Fi + rro)(F2C2 + F3C3) + 2F2F3^, (D4) 



where 



X ■ 
y ^ 



l:!_l+l±J+g_Q^ 

1 + A 2ijji^i^ 

(72 ,72 ? + 73 fla q'MK.ql^q' , 
(72 ,72 9' - 73 ,73 q')/l^qK^' , 
(72 ,73 q' + 73 ,72 g')/KqK^ , 



(D5) 



and 



C2 



C3 



72 g- 
l^q 
. 73£ 



l2c 



+ 



(D6) 



Intcrband orbiton transitions arc represented by a sec- 
ond term in Eq. (D3). To obtain M and A, one should 
just replace — > —k^ in above equations. (This also 
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FIG. 15. Solid line: The reduction of the magnon energy 
(in units of txJse), calculated from Eq.(D3). Broken line 
is a function 3((5J)(1 — 71 p-) with 63 = 0.2, showing that the 
effect of Scjp can fairly be regarded as an effective reduction 
of NN spin coupling. 



leads to — * W2q» in A and r). One can verify that 

the function Slu^ has a cubic symmetry in a momentum 
space. This property is guaranteed by the high symmetry 
of the underlying orbital order in any level of approxima- 
tions. The longer-range, next-NN spin couplings might, 
of course, be dynamically generated by orbital fluctua- 
tions. 

We show the numerical result for magnon renormaliza- 
tion Swff in Fig. 15. For magnon dispersion in Eq.(D3) 

we used |J| = 0.03{riJg'^) ~ 0.05(riJsf;) as obtained 
from Eq. (73), while orbiton dispersion is calculated with 
H — 1.41 which gives orbiton gap tiJse- As estimated 
in the main text, such a gap would be induced by non- 
diagonal hopping To see a deviation of the magnon re- 
mormalization 5uj from the NN Heisenberg form, broken 
line shows a function 3(JJ)(1 — 71 g-) with 5J = 0.2. It is 
noticed that 5 J is indeed close to dJ^ff = 2/9 obtained 
in the main text for the reduction of NN spin couplings. 
Slight deviations from simple NN model are however vis- 
ible, in particular a stronger softening at the (tttttt) point 
is seen. By a numerical fitting, these deviations can be 
traced back to the appearance of longer-range ferromag- 
netic couplings J2 — —0.003, J3=0, and J4 = —0.007 
(all in units oiriJsE)- This is understood due to longer- 
range orbital singlet correlations along the cubic direc- 
tions, as J4 corresponds to a ferromagnetic coupling be- 
tween second- nearest-neighbor spins along cubic axes. In 
principle, these corrections could be observable as a slight 
enhancement of the ratio of small momentum spin stiff- 
ness to full magnon bandwidth as compare to the NN 
Heisenberg model. We should notice however that more 
quantitative predictions are not possible at the present 
stage of the theory for rather obvious reason: A dynam- 
ical spin-orbital interaction is strong (with the coupling 
constant being of the order of 1), so a more elaborate 
treatment is needed to quantify the strongly correlated 
model under consideration. 
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APPENDIX E: EXAMINATION OF SPIN 
INTERACTIONS IN THE PREVIOUS ORBITAL 
MODELS FOR YTIO3 

In this appendix, we examine the magnetic interac- 
tions in the orbital state previously reported in Refs. 15, 
23- 26, in order to check whether this state can explain 
recent neutron-scattering results on the spin couplings, 
spin canting, and magnon gap. 

The Heisenberg spin exchange coupling: First, we dis- 
cuss bond dependence of the isotropic spin interactions. 
The orbital state reported by Hartree-Fock^^ and band- 
structure^"^ calculations are expressed as 



(El) 



with Uc being an occupation of the xy orbital. Using 
these wave functions, it is easy to obtain from Eq. (2) 
spin exchange couplings along c-, a and &-axes: 

J^''> = ^JseUti + r2r3)(l - n,) - (n - r2)](l - n,) (E2 



and 



(ri -I- r2r3)nl 



-in 

2^ 



^2)(1 



(E3) 



The exchange interactions are presented as functions 
of rj in Fig. 16 for different values of Uc- The "meeting" 
points, where J'^'^^ — J^°'''\ are shown by circles for each 
ric- One finds that the isotropy point for the state with 
He = 0.5 (suggested in Refs. 15, 23) is right at the bor- 
der J] = 0, but the exchange coupling is of the AF sign 
there. For larger Hc one may obtain the isotropy point 
with F coupling, but this requires too large values of rj. 
Moreover, the "meeting" point is extremely sensitive to 
both r] and Uc and requires fine tuning. As shown in the 
inset of Fig. 16, the J'^'^''^ / J*^"^) ratio drastically changes 
even at a small (just within ±5%) variation of Uc, and 
may even reverse the sign. 

Spin anisotropy interactions: Next, we investigate the 
effects of spin anisotropy interactions in the state, 
Eq. (El). We denote the occupied orbital state described 
by wave functions in Eq. (El) by a, while the lowest un- 
occupied state is called (3. The state (3 has wave func- 
tion which is a counterpart of that for a state; e.g., 
iV'i) = VI ~ 'nc\dxy) — ^/rL^\dxz) on site 1. Level sep- 
aration between a and (3 is introduced as A^^. Active 
components of the angular momenta at sites 1, 3 and 2, 4 
are Ix and /j,, respectively. These are expressed in terms 
of orbital doublet operators as follows: 



{lx)i,3^±i{P^a-a^P), 



and 



{ly)2A = TiiP^c 



(E4) 



(E5) 
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FIG. 16. Spin exchange interactions in the model (El) as 
a function of the Hund's coupling rj for different values of the 
a;2/-orbital occupation Uc. Energy unit is Jse = 4:t^ /U. Filled 
circles show the position where spin couplings are isotropic for 
fixed value of ric. Inset: The ratio j'"''/ j''^' as a function of 
ric for fixed 77 — 0.24. 



As usual, the spin anisotropy Hamiltonian follows from 
perturbation theory with respect to spin-orbit coupling 
Hso and superexchange interaction Hse- In the SE op- 

"(c) " (c) 

erators J^j and K^j, the following terms contribute to 
the spin anisotropy Hamiltonian: (i) the terms operating 
in the P excited states, and (ii) the terms connecting a 
ground state a with excited ones (3. 

Symmetric spin anisotropy interaction: Third-order 
perturbation with respect to Hgo and Hse gives sym- 
metric spin anisotropy Hamiltonian Hani- For 1-3 bond 
along the c-direction, Hani is given by 



Hani = M^lxSsx), 



where 



A = JsE ^^^/^ fir~) 4nc(l-nc). 
4 VAq^/ 



(E6) 



(E7) 



Here, the terms proportional to small numbers (ri — r2) 
and (r2— ^3) are neglected. For 2-4 bond along c-direction 



H, 



(24) 



A{S2vS. 



2yOiy) 



(E8) 



Note that spin components correspond to that of active 
angular momentum and the interaction is of the AF sign. 
This originates from the fluctuation of active angular mo- 
menta with AF correlation between NN sites, as can be 
inferred from Eqs. (E4) and (E5). For the 1-2 bond in 
the a6-plane, Hani has the following form: 



H. 



(12) 



1 



1 



A{SlyS2 



(E9) 



The combination of spin components in the last term 
is different from that of active momenta at sites 1 and 
2. Such terms originate from the following processes 
(combination of operators): {lixSix{SizS2z)hyS2y)orb, 
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FIG. 17. Antisymmetric DM interaction pattern obtained 

for the orbital model (El). Gray arrow denoted by dij shows 
the orientation of DM vectors. Black arrows represent the 
direction of spins favorable for these interactions. 

{hyS2y{SizS2z)hxSix) orb-, stc. Fov 3-4 bond, symmet- 
ric spin anisotropy Hamiltonian Hani has the same form 
as Hj^J^), where site index 1 (2) is replaced with 3 (4). It 
can be shown that Hani does not cause finite spin canting 
when I J| > A/ 2, where J is the isotropic spin couphng. 
However, it leads to a large magnon gap of classical ori- 
gin, which wc obtained to be ^/SAS. 

DM interaction: For the 1-3 bond, second-order per- 
turbation with respect to Hg^ and Hse gives 

^DM = ^^'^^^13 • (^1 X S3) , (ElO) 

where constant D^'^^ is obtained as 

5^^^ = Jse'^^K^HI - n^f'. (Ell) 

DM vector di3 is given as ^13 = (—1,0,0). We point 
out here that DM constant D'^"^ is much larger than 
D in Eq. (87), obtained in the main text for the SE- 

drivcn orbital states. The reason is that the orbital or- 
der given by Eq.(El) has not that high symmetry, and 

the terms proportional to (ri + r^) in the operator 
do contribute to DM interaction. Therefore, the ratio 
D'^'^^/D DC ocl/j? is large. 
On the in-plane bond 1-2, we find 

H^DM = lD(^'^di2 ■ (5i X §2) , (E12) 

with = 5(^)nc/(l-nc) and di2 = (-1, -1,0). The 
DM interactions on the 2-4 and 3-4 bonds arc given by the 
same forms as Eqs. (ElO) and (E12), respectively, where 

d24 = (0,1,0) and ^34 = —di2- The DM interaction 
and related spin structure in the orbital model (El) are 
schematically shown in Fig. 17. 

Spin caiiLnuj and niucjn.on (jap: First we estimate 
anisotropy constants. Wc consider the orbital state with 
Uc = 0.6 which has a chance to explain the J*^"^-* = J^"*") 
property, giving isotropic spin coupling J about —O.IJse 



at 77 = 0.24 (sec Fig. 17). Given these parameters, one 
obtains from Eqs. (E7) and (Ell), 

I~11(A/A„^)2|J|, 
~ 1O(A/A„0)| J|, 
D'-''^) = (3/2)5("). (E13) 

By minimizing classical energy of DM and isotropic 
Heisenberg interactions, we find that spins cant away 

from the c axis by an angle d ~ 3.5(A/A„/3). 

Also, within a linear spin-wave theory we estimate a 
magnon gap generated by symmetric anisotropy and DM 
interactions as VSAS ~ 9.5(A/A„/3)2| J| and 4S\J\e'^ ~ 
25{X/ Aafj)"^],!], correspondingly. With A ~ 19 meV and 
I J| ~ 2.75 mcV, wc obtain that the ratio (A/A^,^) must 
be less than 0.05-0.06 in order to be consistent with the 
observed canting angle (~ 0.17 rad) and an upper limit 
for the gap (~ 0.3 meV). Thus, the splitting of the lowest 
orbital doublet Aq/j should be at least about 300 meV. 
However, this is hard to reconcile with almost equal four 
short Ti-0 bonds in titanates suggesting an almost de- 
generate doublet picture. 

Based on the above analysis, we think that the orbital 
state (El) predicted by band-structure calculations is not 
supported by recent neutron-scattering experiments in 
YTiOs.^" This is perhaps not really surprising, as a Mott 
insulator with orbital degeneracy represents a strongly 
correlated system, which is difficult to address in a frame- 
work of weakly interacting electrons. 
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